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1.  STATEMENT  OF  PROBLEM 


Th«  research  supported  by  this  contract  was  concerned  primarily  with  (1) 
the  electromagnetic  modes  and  (2)  the  radiation  pattern  characteristics  of 
contemporary  injection  laser  structures.  The  objective  of  this  work  was  to 
investigate  laser  structures  which  can  produce  stable-beam  radiation  patterns. 
Because  the  radiation  pattern  in  the  plane  of  the  p/n  junction  of  many 
injection  devices  is  Influenced  by  the  cavity  geometry,  stable-beam  patterns 
can  be  obtained  only  fr on  devices  with  proper  and  efficient  control  of  the  resonant 
cavity  modes. 

The  control  of  lateral  fields  in  lasers  has  been  eccompllshed  with  many 
novel  devices  such  as  the  channeled-substrate-planar ,  burled  heterostructure 
and  strip-loaded  waveguide.  These  devices  tend  to  be  extremely  complicated 
technologically  so  that  it  Is  advantageous  to  produce  structures  such  as  the 
simple  stripe  geometry  device.  The  present  drawback  to  the  stripe  laser  is  that 
it  is  laterally  unstable.  Much  of  our  work  has  been  aimed  at  rigorously 
modeling  the  electromagnetic  modes  of  stripe  lasers.  We  have  also  applied  our 
analysis  to  other  lasers  such  as  the  channeled-substrate~planar  device. 


2.  SUMMARY  OF  RESULTS 


The  rapid  development  of  semiconductor  laser  technology  since  the  mid  1960's 
has  pushed  the  field  of  optical  communications  to  high  levels  of  promise.  Early 
modeling  of  the  electromagnetic  modes  character  in  the  plane  perpendicular 
to  the  p/n  junction  was  successful  because  the  optical  fields  were  confined 
due  to  the  incorporation  of  optical  barriers  at  heterojunctions. ^  Original 
cavity  thicknesses  distance  between  heterojunction#  were  d  -  1  um  (Fig.  1). 


Fig.  1  Cross  section  of  a 
stripe  geometry  laser  structure. 
The  active  layer  is  bounded  by 
the  two  AlGaAs  layers. 


In  an  effort  to  reduce  the  current  density  at  threshold!  the  cavity  thickness 
was  reduced  to  values  d  <  0.2 w  m.  These  small  thicknesses  play  a  major  role 
in  the  Influence  of  the  lateral  mode  behavior.  The  intense  electric  field  in 
the  cavity  tends  to  stimulate  recombination  of  electron/hole  pairs  and  thua 
depletes  the  carrier  density  at  the  peak  points  of  the  optical  field.  To  put 
this  phenomena  in  perspective  we  must  understand  the  mechanisms  producing 
lateral  confinement  of  the  optical  fields.  Optical  mode  guiding  along  the 
lateral  direction  occurs  because  of  changes  in  the  real  and  Imaginary  parts 
of  the  dielectric  constant.  In  stripe  lasers,  the  optical  mode  is  confined 
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to  regions  directly  below  the  stripe.  Thus  the  dielectric  constant  below  the 
stripe  differs  from  its  value  in  lateral  adjacent  regions  bordering  the  stripe 
region.  For  example,  in  Fig.  2  we  show  a  cross  section  of  a  specific  laser; 


Fig.  2  Stripe  contact  laser  with  active  region  width  V  . 

Mode  gain  occurs  in  the  active  region  while  losses  take 
place  in  the  passive  regions. 

the  active  region  is  below  the  stripe  while  the  passive  regions  are  below  the 
oxide  layer.  Factors  affecting  the  dielectric  distribution  are: 

1)  Cain  in  the  active  and  loss  in  the  passive  regions, 

2)  Lateral  temperature  variations, 

3)  Carrier  density  variations. 

When  the  peak  optical  field  depletes  carriers  under  the  stripe,  the  gain 
changes  accordingly.  Thus  the  refractive  index  changes  according  to  (1)  and 
(2)  above. 

In  our  work  we  have  shown  how  asymmetric  modes  can  be  excited  when 
the  dielectric  distribution  becomes  distorted. ^  These  asymmetric  modes 
cause  extreme  distortions  in  the  radiation  pattern.  Furthermore,  the  amount  of 
distortion  Is  current  dependent  which  reflects  perturbation  in  the  dielectric 
profile. 


Another  major  accomplishment  of  our  work  is  the  first  rigorous  modeling 

(3  4) 

of  the  cavity  modes.  '  All  previous  analyses  of  the  electromagnetic  modes 
of  stripe  lasers  have  been  approximate.  Ue  have  shown  that  the  vertical  and 
lateral  modes  of  cavity  structure  are  interconnected  and  that  a  complete 
boundary  value  solution  of  the  lateral  modes  must  include  interaction  wi th 
vertical  field  distributions. 

Uith  laser  fabrication  technology  developed  at  lt6  present  level,  it  is 
possible  to  reproducibly  grow  thin  epitaxial  layers  such  as  those  in  Fig.  1. 
Consequently,  it  was  necessary  to  explore  the  modal  effects  of  the  lasing 
field  extentions  into  the  various  layers  of  Fig.  1.  In  our  work  we  have 
accomplished  an  extensive  analysis  of  the  effects  of  device  geometry  on  lateral 
mode  behavior  in  stripe  contact  and  channeled-substrate-planar  lasers. 

As  mentioned  previously,  the  problems  of  lateral  mode  instabilities 
in  stripe  contact  lasers  have  been  paramount  because  there  i6  a  lack  of 
lateral  index  guiding.  Optical  stripline  structures  have  a  "built-in"  lateral 

index  step/  '  Consequently,  we  have  modeled  the  stripline  in  an  effort 

(8) 

to  evaluate  Its  mode  behavior  as  an  injection  laser. 

Finally,  we  have  devised  an  experimental  technique  for  estimating 

(9) 

refractive  index  steps  of  passive  waveguides.  The  intent  of  this  effort 

was  aimed  at  evaluating  mode  behavior  of  a  laser  structure  in  its  wafer  form. 
After  the  laser  wafer  is  processed,  slivers  are  cleaved  (the  width  of  the  cleave 
actually  forms  the  laser  length).  The  slivers  are  then  sawed  to  form  the 
actual  devices  for  mounting.  Using  our  technique  the  waveguldlng  properties 
of  the  lasers  in  sliver  form  can  be  evaluated  before  final  processing. 
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Introduction 


The  research  carried  out  by  this  project  is  aimed  at 
understanding  the  tvo~diaensional  radiation  characteristics  of 
semiconductor  injection  lasers.  Basically,  the  radiation  pattern 
characteristics  of  injection  lasers  can  be  divided  into  tvo 
categories:  (1)  the  pattern  in  a  plane  perpendicular  to  the 

Junction  plane  and  (2)  the  pattern  in  the  plane  parallel  to  the 
jwctlon.  Much  of  the  previous  work  in  the  literature  is  concerned 
with  the  pattern  characteristics  in  the  plane  perpendicular  to  the 
junction.  This  is,  of  course,  related  to  the  fact  that  the  mode 
structure  perpendicular  to  the  Junction  (the  x-direction  in  Fig.  1) 

e 

can  be  easily  modeled  using  multilayer  waveguides.  On  the  other 
hand,  very  little  modeling  of  the  transverse  V  modes  (y-dlrectlon  of 
Fig.  1)  has  been  made.  This  is  primarily  due  to  the  fact  that  the 
optical  field  confinement  along  the  y-direction  is  not  well  under¬ 
stood  for  some  laser  devices.  In  addition,  the  W  mode  selection 
characteristics  are  not  well  understood.  For  example,  some  stripe 
geometry  lasers  operate  in  high-order  transverse  W  modes,  yet  the 
high-order  modes  have  less  optical  feedback  than  some  of  the  low  orders. 

In  an  effort  to  understand  the  waveguide  characteristics  of 
the  transverse  V-modes,  we  are  modeling  the  radiation  pattern  in  the 


Research  Direction 

The  simple  Fourier  transform  formula  to  connect  the  far- 
fleld  radiation  pattern  and  the  near  field  haa  been  used  to  analyze 
the  transverse  direction  radiation  phenomena.  In  moat  cases,  the 
aeasured  pattern  is  soaevhat  sharper  than  the  theoretical  one  and  at 
0  ■  ♦  90*  (along  the  x-axis  of  Fig.  1),  the  siaple  Fourier  transform 
formula  does  not  lead  to  a  null  point  at  the  facet.  Froa  ray  optics, 
Snell's  lav  indicates  there  will  be  a  null  at  +  90*  due  to  the 
total  internal  reflection.  Butler  et  al . ^  represented  a  aode  by  a 
spectrum  of  plane  waves  and  determined  the  transmission  characteristic 
of  each  plane  wave.  They  compared  calculated  values  with  experimental 
patterns  to  determine  the  device  geometrv  and  dielectric  parameters 
of  various  layers.  Hockham6  used  a  rather  complicated  procedure  to 
obtain  a  formula  which  added  a  new  factor  to  the  Fourier  transform 
formula.  Levin  7  interpreted  this  term  as  "Huygens  obliquity  factor". 
Both  Butler's  and  Hockham's  formula  applied  the  saddle-point  asymptotic 
technique  to  evaluate  integrals  and  a  cos  0  term  was  obtained  to 

g 

satisfy  the  Snell's  law  at  0*  +  90*.  Recently,  Levin  Included  a 

a  multiplier  into  the  radiation  equation  by  using  a  weighted  average 

resulting  from  the  mode- conversion  effects  occurring  on  the  laser-air 

Interface.  This  method  is  relatively  accurate  and  also  can  be  used  to 

9 

calculate  modal  reflectivity  which  compares  well  with  Ikegami's 
nimterlcal  calculation.  Reinhart  et  al. ^  proposed  a  method  for  cal¬ 
culating  the  modal  reflectivity  by  decomposing  the  modal  fields  into 
an  infinite  masher  of  plane  waves;  the  plane  wave  reflection  coefficient 
was  used  to  obtain  the  modal  reflection  coefficient. 


The  above  considerations  are  baaed  on  a  one  dimensional  model. 


It  la  assigned  that  the  facet  fields  are  uniform  along  the  junction 
(lateral  direction).  The  research  in  this  project  is  concerned  with 
the  model  structure  in  the  lateral  direction  which  relates  also  to 
the  radiation  pattern  in  the  plane  parallel  to  the  junction.  Con¬ 
sequently  ,  a  model  of  the  laser  cavity  is  being  developed  to  model 
the  lateral  cavity  dimensions.  This  model  will  be  applicable  to 
presently  fabricated  devices  such  as  stripe-geometry  lasers11,  buried 
heterostructure  lasers12,  and  Internally  striped  planar  Users.13 

Gordon extended  Reinhart's1^  model  from  one  to  two  dimensions 
which  covers  the  finite  waveguide  width  in  junction  plane,  however,  no 
explicit  relation  between  the  reflectivity  and  the  stripe  width  has 
been  explored.  Naah1^  assumed  that  the  gain  under  the  stripe  contract 
and  the  loss  outside  the  stripe  will  produce  a  complex  dielectric 
constant  and  gain  with  parabolic  shape  along  the  junction  plane  due  to 
the  current  spreading  under  the  stripe.  Experimental  evidence  was  reported 
by  Cook  and  Nash1^  for  gain-induced  guiding  and  refractive- index 
antiguiding  along  the  junction  plane  of  the  stripe-geometry  DH  lasers. 

In  the  direction  normal  to  the  junction  plane  (transverse  direction), 
the  fields  are  confined  by  the  index  steps  between  the  active  region  and 
the  surrounding  passive  medium  and  have  the  sinuslodal  shape.  The  modes 
can  be  restricted  to  the  lowest  order  by  staking  the  active  region  very 
thin  or  index  step  very  small.  However,  along  the  junction  plane,  the 
SK>de  structure  remains  somewhat  uncontrollable  despite  the  introduction 
of  the  stripe  contact.  We  have  investigated  the  two  dimensional  radiation 
pattern  and  mode  reflectivity.  The  active  region  under  the  stripe  is 


treated  as  an  antenna  aperture  and  a  rectangular  waveguide  aodel  la 
applied  to  atudy  the  radiation  and  nodal  reflectivity.  Although, 
the  change  of  index  along  the  lateral  direction  is  continuously 
varying  instead  of  abruptly  changing,  the  gradient  of  the  refractive 
index  can  be  neglected  when  the  index  variation  is  slow.  The 
transaltted  field  can  be  obtained  by  matching  appropriate  fields  at 
the  laser-air  Interface.  The  half-power  beasvldth  in  the  transverse 
direction  conpares  well  with  experiments  while  the  lateral  half-power 
beamwldth  has  a  narrower  value  than  the  measured  one.  This  may  be 
due  to  our  neglecting  the  loss  outside  the  stripe  in  the  active  region. 

The  two  dimensional  modal  reflectivity  has  been  obtained  by 
weighted  average  of  the  electric  field  with  a  specific  multiplier  to 
include  the  effect  of  mode  conversion.  The  relation  between  the 
reflectivity  and  the  stripe  width  is  shown.  The  fundamental  lateral 
mode  has  the  largest  re f lectivi ty  for  typical  stripe  widths  except  near 
the  cut-off  width  of  the  higher  lateral  modes.  For  large  stripe  widths, 
all  modes  approach  the  asymptotic  value  which  can  be  obtained  from  the 


one  dimensional  model. 
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Pig.  1  Radiation  from  the  end  facet  of  an  Injection  laser.  The  facet  Is  located 
In  the  x-y  plane  with  the  2-axis  normal  to  the  facet.  The  half-power 
beam  widths  6xand  are  measured  In  the  planes  perpendicular  and 
parallel  to  the  junction  plane. 
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Introduction 


The  research  performed  on  this  project  has  been  directed  toward 
the  analysis  of  semiconductor  laser  modes  in  the  plane  of  the  Junction. 

(Mode  confinement  perpendicular  to  the  active  region  is  due  to  the 
grown  heterojunctiom  while  confinement  along  the  junction  plane  is 
associated  with  the  device  geometry.  In  the  case  of  the  stripe-geometry 
laser,  confinement  is  due  to  the  stripe  contact  as  shown  in  Fig.  1.) 

More  specifically,  we  give  brief  results  of  our  work  concerned 

with:  (1)  a  new  waveguide  model  developed  for  stripe  geometry  lasers, 

(2)  a  comparison  of  theoretical  and  experimental  radiation  patterns  in 
the  plane  of  the  Junction  (this  pattern  is  tied  to  the  lateral  waveguide 
geometry),  and  (3)  two  dimensional  mode  reflectivity  data. 

Research  Direction 

In  Fig.  1(a)  we  show  a  cross  section  of  the  laser  geometry  while 

in  (b)  we  show  the  detailed  geometry  of  the  active  region.  In  the 

direction  normal  to  the  Junction  plane  (transverse  direction),  the  fields 
are  confined  by  the  index  steps  between  the  active  region  and  the  sur¬ 
rounding  passive  media  and  have  the  sinusoidal  shape.  The  mode  car.  be 
restricted  to  the  lowest  order  by  making  the  active  region  very  thin  or 
the  index  step  very  small.  Hcwever,  along  the  Junction  plane,  the  mode 
structure  remains  somewhat  uncontrollable  despite  the  introduction  of 
the  stripe  contact.  The  active  region  under  the  stripe  is  treated  as  an 
antenna  aperture  and  a  rectangular  waveguide  model  is  applied  to  study 
the  radiation  pattern  and  mode  reflectivity. 

Figure  2  shows  the  far-field  radiation  pattern  in  both  the  transverse 
(vertical)  and  lateral  (horizontal)  directions.  The  structure  has 
n,  •  3.6,  n^  •  n^  •  3.46,  n^  ■  n^  "  3.595,  d  •  0.25  um,  and  2w  •  13  urn. 

The  vertical  direction  radiation  pattern  is  shown  in  Fig.  2  (a)  where 
the  half  power  beamwidth  is  40.5*  and  is  roughly  the  same  as  the  one 
obtained  from  one  dimensional  model .  The  lateral  direction 
radiation  pattern  has  *  3.6*.  This  is  somewhat  narrow  than  the  measured 

result  which  is  shown  in  Fig.  2  (b) .  This  discrepancy  may  be  due  to  the 
uncertain  estimation  of  n^  and  n^  outside  the  active  region  in  lateral 
direction.  Also,  the  difference  may  be  due  to  the  fact  that  losses  in 
regions  4  and  5,  outside  the  active  region,  were  neglected. 


To  more  accurately  model  the  radiation  pattern  in  the  horizontal 
direction,  a  waveguide  geometry  different  from  that  shown  in  Fig.  1  must 
be  explored.  Current  spreading  under  the  stripe  contact  primarily  deter¬ 
mines  the  field  distribution  in  lateral  direction.  A  descrepancy  between 
measured  and  calculated  results  may  be  due  to  the  fact  that  current 
spreading  was  neglected.  The  refractive  index  does  not  change  abruptly. 
Recently,  Cook  and  Nash*  claimed  they  found  experimental  evidence  to 
show  that  lateral  field  confinement  is  due  to:  (1)  gain-guiding,  and 
(2)  index  anti-guiding.  Overall,  the  gain  Induced  gu  once  is  predominant. 

In  order  to  approximate  the  gradual  change  of  the  index,  we  apply  the 
2 

five-layer  model  to  represent  the  lateral  direction  change.  For 
simplicity,  only  three-layer  structure  is  considered  in  vertical  direction. 
The  model  of  the  waveguide  structure  is  shown  in  Fig.  3.  We  as sun?  that 
the  center  region  3  has  strong  gains  50  cm  1  with  large  thickness; 
regions  2  and  4  have  less  gains  20  cm  *  with  narrow  widths.  However,  the 
outer  regions  1  and  5,  outside  the  stripe  contact,  have  losses  -10  cm 
In  order  to  Include  the  anti-guiding  effect,  we  assume  that  the  center 
layer  has  the  lowest  index  3.599  in  lateral  direction;  two  surrounding 
regions  2  and  4  have  indices  3.5995;  and  two  outer  regions  1  and  5  have 
indices  3.6.  The  vertical  confinement  comes  from  the  heterojunction  index 
steps  in  which  we  have  n^  •  a_  •  3.46,  n^  ■  3.599. 

By  suitable  width  adjustment  of  various  regions  and  choosing  gain- 
loss  values,  we  can  fit  the  calculated  lateral  radiation  pattern  to  a 
measured  one  very  accurately.  For  example,  we  choose  d^  ■  d^  ■  2.625  um, 
d,  -  8.75  um  and  the  index  and  gain/ loss  values  of  various  regions  as 
above  in  which  the  gain  guiding  and  index  anti-guiding  effects  are  con¬ 
sidered.  The  result  is  shown  in  Fig.  4.  The  calculated  half  power  beam- 
width  is  7.2*  which  is  about  the  same  as  the  measured  one. 

Finally,  the  modal  reflectivity  of  various  two-dimensional  modes 
has  been  obtained  as  a  function  of  the  stripe  width.  The  results  are  shown 
in  Fig.  5.  (The  parameters  are  listed  in  the  figure.)  We  notice  that  the 
fundamental  lateral  mode  has  the  largest  reflectivity  except  near  the  cut¬ 
off  width  of  the  higher  lateral  modes.  This  means  that  the  lower  order 
lateral  modes  have  a  better  chance  to  oscillate  than  the  higher  order  ones. 
Also,  there  are  minimal  values  near  the  cut-off  stripe  width  of  the  various 
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lateral  modes.  As  the  stripe  width  is  reduced  further,  the  reflectivity 
of  higher  orders  increases  very  sharply  and  may  have  a  value  larger 
chan  the  fundamental  mode.  Generally,  the  lateral  fundamental  mode  has 
the  highest  reflectivity.  On  the  other  hand,  at  large  stripe  widths 
R^  R^,  the  one  dimensional  reflectivity. 

References 

1.  D.  D.  Cook  and  F.  R.  Nash,  J.  Appl.  Phvs.  1660  (1975). 
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-  ACTUAL  INDEX  DISTRIBUTION 

APPROXIMATE  F!Ve  LAYER  MODEL 


Five  lave;  model  is  used  to  approximate  Che  lateral 
direction  field  distribution.  The  gain-guiding  and 
lndex-antlguidlng  are  shovn. 
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Fig.  4.  The  far  field  radiation  pattern  obtained  froa  the  five  layer  andel 
calculated  and  aeaaured  reunite  are  in  good  agreeaient. 
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BRIEF  OUTLINE  OF  RESEARCH  FINDINCS 


Th*  work  during  chi«  previous  period  he*  been  aimed  at  developing  a 
comprehensive  model  of  the  etrlpe  geometry  laaer  structure.  Typical 
laaer  wafer*  are  fabricated  by  growing  multiple  layer*  of  (Al  Ca)A* 
with  varying  percentages  of  Al  and  a  final  layer  of  CaAa.  Contemporary 
structure*  have  very  narrow  active  layers  0.1  -  0.2  um)  with  small 
amounts  of  Al  (*v«  TUT).  The  active  layer  is  generally  sandwiched 
between  two  (Al  Ca)A*  layers  containing  approximately  ’ll  Al.  The 
final  LPE  growth  layer  is  CaA*  which  is  used  to  facilitate  the  forming 
of  metallic  contacts.  However, first  S102  i*  deposited  over  the  grown  GaAs 
layer.  Narrow  stripes  of  Sl(>2  are  removed  to  form  th*  light  propagation 
axis.  Finally,  the  metal  layer  is  deposited  over  the  SiC^  and  CaA a  stripes. 

The  stripe  metallic  contact  to  the  CaA*  cap  thus  forms  a  path  of  current 
flow  through  the  grown  layers.  The  object  of  th*  stripe  is  to  limit  th* 
lateral  current  spreading  and  to  laterally  confine  the  recombination  region 
in  the  thin  active  layer. 

The  mathematical  model  which  we  ate  presently  developing  will  allow  for 
a  complete  analysis  of  mode  characteristic*  in  stripe  lasers.  This  model 
takes  in  to  account  the  complete  details  of  the  active  region,  the  grown 
confining  layers,  the  CaA*  cap,  and  finally  the  SiO,  -  metallic  layer 
formed  over  the  GaAs.  We  have  developed  design  curves  which  correlate 
the  number  of  lateral  waveguide  modes  with  the  stripe  width  and  th*  level 
of  pumping  of  the  active  region.  The  calculations  show  that  lateral 
mode  stability  becomes  a  problec  when  lasers  must  be  pumped  at  high  levels. 
Lasers  designed  for  single  mode  operation  must  be  putted  at  low  levels. 

The  complete  details  of  our  analysis  are  in  preparation. 
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BRIEF  OUTLINE  OF  RESEARCH  FINDINCS 


Much  of  the  work  over  the  previous  period  has  been  concerned  with  a  study  of 
modes  in  stripe  geometry  injection  lasers.  It  has  been  found  experimentally  that 
oxide  stripe  lasers  frequently  produce  radiation  patterns  with  major  lobes  cocked 
at  slight  angles  (0-10*)  with  respect  to  the  normal  of  the  end  facet.  In  addition, 
the  lobe  direction  is  a  function  of  drive  current.  (Obvious  application 
considerations  involve  coupling  laser  output  power  to  optical  fibers.  It  is 
Important  to  reduce  beam  Jitter  to  obtain  optimum  coupling.)  In  an  effort  to 
understand  the  lateral  patterns  we  have  modelled  stripe  lasers  with  ar.  asymmetric 
dielectric  profile.  The  theoretical  patterns  of  our  model  also  produce 
asvmmetric  patterns.  Consequently,  we  conclude  that  pattern  instabilities 
are  produced  bv  laterally  changing  dielectric  profiles. 

Finally,  most  of  the  electromagnetic  models  of  injection  lasers  with 
stripe  contacts  are  insufficient  in  their  descriptions  of  lateral  modes. 

We  have  developed  a  new  mathematical  model  useful  for  analysing  the  field 
structure.  For  example,  assuming  a  lateral  parabolic  dielectric  profile 
the  waveguide  modes  are  written  as  linear  combinations  of  He rmlte-Caussian 
functions.  Previous  analyses  have  assumed  that  lateral  modes  are  described 
by  single  Herralte-Caussian  functions.  For  a  giver. parabolic  profile,  the 
new  model  shows  that  the  lateral  field  distribution  of  the  fundamental  mode 
spreads  more  than  that  obtained  using  previous  models.  Consequently,  the 
radiation  beanwidth  will  decrease.  With  the  more  accurate  description  of 
the  lateral  modes  we  can  then  obtain  accurate  dielectric  profiles  which  arc 
governed  by  the  device  geometry. 
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ASYMMETRIC  MODES  IK  OXIDE  STRIPE  HETEP.OJUKCTIOK  LASERS 


J.  K.  butler  and  H.  S.  Scanner s ,  Jr. 

RCA  laboratories 
Princeton,  K.  J.  08540 

ABSTRACT 

Tne  region  between  the  pair  ol  he tero junctions  of  an  oxide-stripe  laser  is 
modelled  by  a  step  profile  in  the  lateral  distribution  ^parallel  to  the  hetero- 
junctions)  of  the  gain  and  of  the  refractive  index,  in  an  attempt  to  pain  insight 
about  the  interpretation  of  the  asymmetric  radiation  patterns  frequently  observed. 

A  region  of  width  V  concentric  with  the  stripe  contact  is  subdivided  into  three 
sections,  in  each  of  which  the  gain  and  index  are  uniform.  Tne  step  heights  for  the 
gain  distribution  and  for  the  index  distribution  are  given  complementary  symmetry, 
with  the  gain  maximum  super  imposed  on  the  index  minimum.  Asvrsaetrv  is  introduced 
bv  lateral  shift  of  the  extrema  at  fixed  V  .  With  values  of  pain,  index  step, 
and  W#  that  seer  reasonable  for  the  laser  studied,  the  neaT  and  far  fields  repro¬ 
duce  the  measured  profiles  of  a  representative  10  um  stripe  laser.  The  offset 
of  the  spontaneous  profile  is  about  3  ue.  f rotr  the  center  of  the  stripe,  and  the 
beam  has  a  single  lobe  for  all  low  order  modes  with  offset  of  order  5'  from  the 
facet  normal.  Because  of  the  near  cancellation  of  gain  guiding  by  index  anti¬ 
guiding  which  is  required  to  give  the  observed  asymmetries ,  the  propagation  con¬ 
stant  becomes  very  sensitive  to  injection  level;  this  may  explain  the  observed 
non-»ono tonic  change  of  modal  power  and  the  very  rapid  shift  of  modal  wave¬ 
length  with  current. 

•Work  supported  in  part  bv  the  U.S.  Arrv  Research  Office. 

♦Permanent  address.  Southern  Methodist  University,  Dallas,  TX  75275 


ABSTRACT 


LATERAL  MODES  OF  STF.IPE  GEOMETRY  INJECTION  LASERS* 

bv 

J.  K.  Butler  and  J.  B.  Delaney 
Southern  Methodlat  I'nlveraltv 
Dallas.  TX  75275 

A  new  eathesat ical  s>nael  useful  for  analvring  lateral  Bodes  of 
stripe  geometo  lasers  is  presented.  The  oxide  stripe  laser  is  aodeled 
as  a  three  iaver  waveguide  in  which  the  dielectric  constant  of  the  active 
I aver  varies  only  along  the  lateral  direction;  the  dielectric  constant 
of  the  surrounding  passive  lavers  is  assumed  ti  be  position  independent. 

The  solution  technique  affords  a  rigorous  Batching  of  the  fields  of  the 
active  laver  with  those  of  the  surrounding  passive  lavers.  To  illustrate 
the  oodel,  the  nodes  of  a  waveguide  with  parabolic  dielectric  variation 
along  the  lateral  direction  are  investigated.  The  fields  are  written  as  a 
linear  coafcination  of  Henrite-Caussian  (K-G)  functions;  heretofore, 
fields  have  been  described  with  a  single  H-C  function.  The  lateral  field 
distribution  of  the  fundamental  mode  spreads  more  than  the  one  obtained 
using  single  H-G  functions.  Fundamental  aode  spread  (spot  sire  at 
halfpower)  is  calculated  and  related  to  the  gain  distribution.  In  addition, 
the  peak  gain  fields  are  determined  at  threshold  for  various  waveguide 
geometries . 


•Supported  bv  the  V .  S.  An-  Research  Office. 
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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


Over  this  last  period  we  have  developed  for  the  flret  time  a  mathematical* 
model  for  analyzing  the  optical  fields  of  stripe  geometry  injection  lasers.. 

The  aim  here  has  been  to  develop  such  a  model  so  that  an  accurate  description 
of  the  lateral  dielectric  profile  can  be  used  for  analyzing  the  lateral 
modes.  Previous  work  on  this  grant  has  shown  the  important  result:  asym¬ 
metric  dielectric  profiles  lasyvetrlc  about  the  center  of  the  stripe  contact) 
lead  to  a  distortion  of  the  mode  radiation  pattern.  With  our  new  rigorous 
model  complete  radiation  pattern  behavior  can  be  understood. 

lr.  this  previous  period  we  have  developed  a  model  which  describes  the  number 
of  lateral  modes  which  mav  propagate  in  stripe  geometry  lasers  and  related 
devices.  This  study  also  addresses  the  problems  of  "kinks"  or  the  high  non- 
linearities  in  the  P-1  diode  characteristics.  These  kinks  are  primarily 
caused  bv  the  drive  sensitive  dielectric  profile  which  of  course  governs  the 
mode  pattern  and  the  mode  content. 

In  conclusion,  we  have  greatly  increased  our  understanding  of  the  performance 
characteristics  of  various  laser  devices.  It  is  clear  that  the  classical 
strlpe-geomecrv  devices  will  have  to  a  certain  degree  mode  stability  and 
pattern  distortion  problems. 


ABSTRACT 


LATERAL  MODI  CONTENT  Of  STRIPE  GEOMETRY  LASER  STRUCTURE?* 

by 

J.B.  Delaney  and  J.K.  butler 
Southern  Mathodlat  University 
Dallas.  Texas  75275 

and 

H.  kressel 
RCA  Laboratories 
Princeton,  N.J.  0854C 

A  model  Is  presented  describing  the  lateral  node  content  cf 
stripe  geooetrv  laser  structures.  A  dielectric  dlscoctmultv  ir 
the  lateral  junction  plane,  determined  froc  the  gain/loss  profile 
of  the  active  laver  and  fros  the  total  geometrical  structure  cf 
stripe  geometrv  devices. Is  calculated  and  related  to  the  cutoff 
conditions  of  the  various  lateral  modes.  In  our  model,  character¬ 
ized  by  a  step  In  the  lateral  complex  dielectric  constant  of  the 
active  laver.  the  cutoff  conditions  are  dependent  upon  the  pair, 
region  width  defined  bv  Injection  under  the  stripe.  The  model 
also  provides  a  comparative  analysis  of  the  standard  oxide  stripe 
and  the  channeled -substrate  lasers.  We  examine  conditions  under 
which  the  fundamental  mode  Is  unstable  and  In  fact  may  not  propagate. 


*Thls  work  was  supported  bv  u.S.  Arwv  Research  Office. 
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A  RIGOROUS  BOUNDARY  VALUE  SOLUTION  FOR  TUI 
LATERAL  MODES  OF  STRIPE  CEOMETRY  INJECTION  LASERS* 

t 


by 

J.  K.  Butler  and  J.  B.  belanev 
Southern  Methodist  Unlversitv 
Dallas.  TX  75275 

A  new  matneaaticai  mood  useful  dor  anaivzmg  lateral  nodes  of 
stripe  geomerrv  lasers  is  presented.  The  oxiae  stripe  laser  Is  modeler 
as  a  tnree  laver  waveguide  In  which  the  dielectric  constant  of  the  active 
laver  varies  onlv  along  the  lateral  direction;  the  dielectric  constant 
of  the  surrounding  passive  lavers  Is  assumed  tc  be  position  independent. 
The  solution  technloue  affords  a  rigorous  matching  of  the  fields  of  the 
active  lever  with  those  of  the  surrounding  passive  lavers.  Tc  Illustrate 
the  model,  the  modes  of  a  waveguide  with  parabolic  dielectric  variation 
along  the  lateral  direction  are  Investigated.  The  fields  are  written  as  e 
linear  combination  of  Hermit e-Gauss lan  (H-C)  functions;  heretofore, 
fields  have  oecn  described  with  a  single  K-G  function.  Fundamental  mode 
spread  (spot  size  at  half power)  is  calculated  and  related  to  the  pair 
distribution.  ! Previous  estimates  of  the  lateral  field  spread  of  the  fund¬ 
amental  mode  using  a  single  H-C  function  not  rigorously  matched  at  the 
boundaries  can  vleld  spot  sises  as  much  as  301  different  from  results 
calculated  frov  linear  combinations  of  H-C  functions.)  In  addition,  the 
peak  gain  fields  are  determined  at  threshold  for  various  waveguide  geo¬ 
metries. 
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BRIEF  OUTLINE  OF  RESEAkO;  F! NDINGS 


’’he  work  over  this  past  period  has  progressed  along  two  lines:  (1)  Theoretical 
investioation  of  lateral  waveguide  modes  and  (2)  Experimental  investigation  of  wave 
propagation  in  laser  structures. 

J.  B.  Delaney  has  worked  primarily  in  theoretical  aspects  of  lateral  mode 
cnaracteri sties .  One  of  the  primary  research  tonics  in  injection  laser  studies 
concerns  the  development  o'  devices  which  support  staple  oscillating  modes 
InstaPi  1 1  ties  produce  "kinks1'  in  the  power-current  cnaracteristics  and  radiation 
patterns  with  lobe  directions  that  are  sensitive  to  drive.  Stripe  geometry 
lasers  have  instabilities  which  are  attributed  to  the  index  antiguiding  nature 
o *  tne  waveguide  modes.  Index  antiguiding  is  primarily  due  to  earner  injection 
in  the  active  region  and  thus  depresses  the  refractive  inaex  there.  We  have 
shown  that  the  overall  device  geometry  of  stripe  lasers  can  contribute  to 
antiguiding  effects 

H.  W.  Scott  has  been  working  on  the  experimental  investigations.  (He  has 
also  made  a  study  of  radiation  patterns  'rom  optical  stripline  waveguides. 

*nese  waveguides  are  constructed  similar  to  stripe  geometry  lasers.  He  is 
in  the  final  stages  o'  his  manuscript  preparation  "Radiation  from  Optical 
Strip' me  Waveguide."  The  experimental  measurements  being  conducted  concern 
tie  focusing  of  1.15  micron  radiation  from  a  He-Ne  laser  on  the  facet  of  a 
stripe  geometry  laser  and  observing  the  mode  structure  on  the  opposite  facet. 

'he  lasers  were  prepared  at  RC*  Laboratories ,  Princeton,  N.  J.  and  were  mounted 
on  special  headers  so  that  observations  of  front  and  rear  facets  could  be 
made.  Our  basic  objective  for  this  study  is  to  develop  an  understanding 
o'  waveauiding  in  lasers  by  probing  the  waveguide  with  light  at  1.15  microns 
which  is  insensitive  to  the  waveguide  gain  which  guides  the  light  at  the 
0.S3  micron  wavelength. 
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BRIEF  OUTLINE  OF  RESEARCH  FINDINCS 


The  work  over  this  past  period  Has  progressed  along:  (1)  theoretical 
investigations  of  lateral  waveguide  modes  and  (2)  experimental  investigation* 
of  wave  propagation  in  laser  structures.  Our  aim  is  directed  toward  the 
development  of  stable,  single  mode  lasers  important  for  many  system  applications. 

J.  B.  Delaney's  work  in  the  analysis  of  laser  modes  is  proving  to  he 
extremely  valuable  for  ctia rac t er 1 t ing  contemporary  loser  structures.  Hi* 
work  discussed  in  the  last  progress  report  (Details  to  be  published  in 
IEFE  J.  Quan.  Electron.,  Aug.  1978)  concerned  the  development  of  a  wave¬ 
guide  model  for  stripe  lasers  with  proper  attention  focused  on  the  overall 
geometry.  Including  the  metallic  stripe  contact  and  isolating  S10,  layers. 

The  analysis  was  also  applied  to  channcl-aubsrrato-planor  (C SD  lasers  which 
are  relatively  mode  stable.  Our  calculations  indicated  that  CSP  devices 
are  more  mode  stable  than  rhe  stripe  iasvrs.  The  analysis  it  being  further 
developed  in  order  to  explore  its  full  implementation  on  new  devices. 

J.  B.  Delaney's  recent  work  which  i-.  being  prepared  for  publication 
concerns  the  far  field  radiation  of  stripe  lasers  with  small  stripe  widths 
S  (<5  urn).  Present  literature  indicates  tne  lateral  field  beam  (parallel 
to  the  Junction  plane)  decreases  as  the  active  layer  thickness  d  decreases 
(below  0.2  urn).  Our  new  work  shows  that  t he  la tera 1  beam  broadens  as  d 
drops.  Further,  the  required  gain  at  threshold,  glh,  is  developed  in  term* 
of  S  and  d.  Our  rest. Its  show  that  for  a  given  d,  an  optimum  S  cat.  be  deter¬ 
mined;  this  correlates  to  the  determination  of  a  minimum,  value  of  current  at 
threshold. 

Much  of  the  analysis  mentioned  above  require  knowledge  of  the  optical 
properties  (the  dielectric  constat.?)  of  the  device  material  (epitaxial  SiO., 
and  metallic  laversl.  The  epitaxial  and  SiO,  lavers  have  published  optical 
properties,  whereas,  those  of  the  metal  contact  are  relatively  unknown. 
(Contemporary  devices  have  active  regions  near  the  contact  with  optical  fields 
penetrating  the  metal.)  Also,  the  mechanism  producing  lateral  mode  confine¬ 
ment  in  stripe  lasers  is  still  an  unknown.  Hence,  we  have  been  developing 
experimental  methods  to  explore  the  various  unknown  parameters.  M.  Scott 
ha*  continued  to  develop  this  important  area  of  our  program.  His  recent  work 
submitted  for  publication  concerned  the  pattern  measurement  of  a  laaar  cavity 
which  was  excited  with  an  external  source.  From  the  far  field  measurement 
we  estimated  the  properties  of  the  epitaxial  grown  layers. 
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Evaluation  of  Dielectric  Optical 

Waveguides  Proa  Their  Far  Field  Radiation  Patterns* 
by  M.  W.  Scott  and  J.  K.  Butler 
Southern  Methodist  University 
Dallas.  Texas  75205 

ABSTRACT 

An  experimental  technique  for  determining  the  characteristics  of 
dielectric  optical  waveguides  is  presented.  The  far  field  pattern  of 
the  waveguide  Is  measured  and  compared  with  a  theoretical  pattern.  The 
presence  of  spurious  light  lr.  the  far  field  complicates  the  oeasureattnt, 
but  this  difficulty  can  be  minimized  for  suiny  structures.  Application 
to  sealconductor  laser  structures  is  discussed.  Determining  the  wave- 
guiding  properties  of  laser  structures  prior  to  processing  can  result 
in  the  early  detection  of  faults  for  Increased  yield. 
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C'ne  class  of  channel  waveguides  frequently 
proposed  *»  elements  ir.  int e g rat f d  optical 
circuit*  irt  optical  strimimes  •  The  struc- 
ture  cor.tltts  of  a r  underlying  substrate  with 
ar.  index  n  j.  a  slab  of  thickness,  d,  having  a 
rafractjvr  index  r.;  ,  and  ar.  ev*rlav*d  strip* 
cf  width.  2w  ar.d  ar.  index  n>.  .  the  Index 
n;  »  r.j  ,  rj  Optical  confinement  cf  the  field* 
;»  restricted  to  the  slab  in  the  region  belos 
the  stripe  Previous  theoretical  investiga¬ 
tion*!. •  of  the**  waveguide*  nave  u»ed  ar. 
artificial  index  discer.tlr.ust'  to  explain  the 
lateral  optical  confinement  Tn*  approach 
described  her*  it  tc  asiunt  the  waveguide 
aodei  are  ccrresed  of  a  continuum  cf  plane 
wave*  Confinement  tc  the  regisr  below  the 
stripe  it  cbfamei  b'  spatial  interference 
The  propagation  eor.startt  are  calculated  as  a 
function  of  waveguide  geometries  Model  cut¬ 
off  conditions  are  determined  as  a  function  of 
the  slab  thickness  d  and  stripe  width  (2w 
for  a  rar.ge  of  index  st^ps  The  observed  nodes 
reported  bv  Blu»  et  al.3  agree  with  cur  calcu¬ 
lations  x>e  have  measured  the  lateral  field 
profiles  of  several  CaA*  optical  striplines, 
exper  imente  1  •  tree  ret  i cal  comparisons  are  ».ad* 

*»  have  fabricated  several  directional 
couplers  consisting  cf  parallel  optical  strip- 
line  waveguides  .he  coupling  lengths  cf 
these  devices  have  been  measured  fer  different 
waveguide  widths  and  separations.  Using  our 
theoretical  model  we  have  calculated  the 
coupling  lengths,  the  experiments!  and 
theoretical  coupling  lengths  are  ir.  agreement. 
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Asymmetric  Modes  in  Oxide  Stripe  Heterojunction 

Lasers 

JEROME  K  BITLER  memmf  ieee.  and  HENRY  S  SOMMERS  JR 
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I  Intooduction 

THIS  paper  prearnti  i  theoretical  model  of  the  lateral 
model  f  model  in  the  function  plane  I  of  the  narrow  oxide 
itnpe  infection  laiei  which  encompaaaet  teveral  expenmental 
feature!  that  are  not  produced  b>  the  popular  Hermite- 
Gauiuar  ( H-G  i  model  1 1  ]  The  investigation  an  prompted 
b>  our  expenmental  itudiei  of  a  large  number  of  oxide  itnpe 
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luen  prepared  frorr.  s  variety  of  duubie-heuro)uncnon  wafers 
with  imr  -aviuei  lr.  particular,  (he  H-G  model  falls  tc  ac¬ 
count  for  (he  followmi  observed  properties  1 )  In  s  lueabk 
fraction  of  individual  lasers  (he  facet  illumination  i  neat-field) 
and  (he  lateral  tfar  fielJi  profile  are  markedh  asymmetncal  in 
respect  to  the  pin  profile  and  the  normal  to  the  facet,  teapot- 
tmly  (2)  Thu  characteristic  has  also  been  noted  in  other  lab¬ 
oratories  (3)  2)  While  the  spatial  patterns  are  often  quite 

complex,  with  much  resolved  structure  they  show  little  vena¬ 
tion  with  current  or  output  power  [2)  Thu  su pests  that  the 
mode  profile  depends  but  weakly  on  mode  otder 

U.  Theoretical  Model 

The  new  theoretical  model  introduces  an  ad  hoc  dieiectnc 
profile  which  pernio  s  mac  rune  study  to  find  the  tan  of  the 
spatiai  isvmmctnes  ir.  the  real  and  imapnars  parts  of  the  di¬ 
eiectnc  constant  to  produce  the  observed  asymmetries  of 
modai  pansms  For  pedagopcai  reasons,  the  treatment  u  car 
in  the  customers  terms  of  index  anoguidmg  and  pin  guidmg 
associated  with  the  free  came;  dutnbunon  but  does  not  im¬ 
ply  that  the  earner  dutnbution  u  the  compute  explanation  of 
the  dielectric  profile  It  u  found  that  i  near -cancellation  of 
the  two  types  of  guiding  topther  with  a  suitable  asymmetry 
common  to  both  produces  the  kinds  of  patterns  observed 
While  the  adjusted  sues  of  the  spatial  asvtnmetnet  are  found 
to  be  consistent  with  the  present  sketchv  ideas  used  to  con¬ 
nect  the  free  carnet  density  and  dielectric  constant,  the  con¬ 
tribution  of  the  paper  u  not  support  or  nontupport  of  their 
utility .  but  rather  knosviedp  of  both  s  profile  which  can  yitld 
such  asvmmetnc  modes  and  of  its  dispersion  relations  [4] 

The  effect  of  the  concave  refractive  index  and  convex  pm 
profiles  u  to  product  modes  confined  laterally  due  io  the  pin 
The  modes  ate  "pin  guided"  (5)-(']  as  oppoaed  to  "index 
guided  "  Now  fot  s  specific  complex  dielectric  profile  related 
to  the  pm  snd  index  profiles,  modai  complex  propeption 
constants  and  field  profiles  come  from  s  solution  of  MaxsveC'i 
equations  Thus,  each  lateral  mode  propaptts  as  expt^.*  > 
where  the  pin  of  the  rth  mode,  u  determined  from  the 
complex  propeptide)  constant  The  value  of  g,  is  in  s  sense 
prorated  over  the  lateral  pm  profile  for  each  mode  It  should 
also  be  noted  that  field  profiles  of  index -guid'd  modes  are 
rather  inae native  to  pm  profiles,  whereas,  field  profiles  of 
gam  guided  modes  are  tightly  coupled  to  the  gam  profile 

The  electromagnetic  field  equations  of  two-dime naonal 
waveguides  can  be  approximated  as  dac timed  by  Marcatfb  (I] . 
and  the  field  equations  developed  by  Meumtng  only  one  trans¬ 
verse  component  of  the  eiectnc  field  The  croa  section  of  the 
laaer  itructsir*  dtosvn  m  Fig  1  ha*  opacai  fields  concentrated 
m  the  vtanitv  of  the  active  region  with  eiectnc  fields  polarized 
parallel  to  the  acme  lave?  Current  ficrwwg  through  the  acne* 
repon  produces  optical  gam  spread  over  t  width  •»,  To  ap¬ 
proximate  the  affect  of  the  nonuniform  lateral  datnbutson  of 
earner  dertarv  the  active  repon  a  partitioned  into  thtee 
-v  gross  which  ire  supported  on  either  nde  bv  passive  regions 
(Fig.  21  Dated  la)  shows  the  paattvi  repons  1  .  ?  and  the  re 
pons  2.  3.  and  4  of  the  acme  section  detail  (bi  (hoses  the  cor¬ 
responding  gams,  which  represent  regeneration  m  the  acme 
mcnon  and  danpatroi  tn  the  paasree  anngr  and  detail  ic) 
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shows  the  index  of  tefracuon.  which  increases  from  the  mini¬ 
mum  in  the  high  gam  central  region.  The  guiding  by  the  gam 
maximum  a  sufficient  to  overcome  the  defocuang  effect  of 
the  index  minimum 

For  arave  propagation  m  the  z  direction  of  the  form 
exptrwr  -  ft)  where  y  •  0.2  *  it.  the  transverse  dependence 
of  the  field  quantities  becomes 
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where  t,  n  the  free  space  wave  number,  n,  •  120*  Cl  and 
•Iz. ,»)  ts  the  re  la  me  dialactnc  constant  The  traneverse  wav* 
funcoon  »  sans  fist 
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where  ?,s  u  the  transverse  Lspiacian  opcrsiot 
As  a  firs:  approximation  we  mumt  that  the  vertical  (a  I  de¬ 
pendence  of  the  wave  function  is  independent  of  s .  the  lateral 
position  in  the  wavtgusding  layer  This  assumption  is  reason 
able  fot  the  lasers  unde:  consideration  which  have  large  di¬ 
electric  steps  at  the  heterojuncuons  compared  to  the  variation 
in  the  dielectric  conitant  m  the  lateral  direction  Come 
quently .  the  solution  to  the  wave  equation  in  the  active  layei 
can  be  written  as 

H*1  *fvt  (7) 

For  the  structure  of  Fig  1.  the  index  step  between  the  active 
layei  and  the  p-  and  n-  confining  walls  are  identical  with  An  =» 
0  1<3  Only  the  fundamental  vertical  mode  propagates  so  that 
*U>  •  cos  k^x  Substituting  *(_v)  into  (t>i  and  simplifsing 

P) 

— *  *0 »(.v >r  ♦  k'i  > .( 1  r>  *i  ♦  t*  W(_v»  •  0 

dy 

(8t 

where 

«(_V)  ■  n’(_v)  -  1  kl  *  i  ft_v!  at  v!  te 

ts  the  lateral  dependence  of  the  dielectric  constant  tn  the  ac 
tive  lavet.  »„  u  the  dielectric  constant  of  the  p  and  n-walli 
and  T  is  the  vertical  confinement  factor  representing  the  frac 
Don  of  the  wave  mode  the  layer  hetweer.  the  heteroiunctions 
P  a  0  6  for  the  structure  of  Fig  1  |o]  The  function  *<_v) 
found  from  (8)  determines  predominantly  the  nest-field  be 
havsor  along  y  and  the  corresponding  lateral  radiauon  pattern 
The  modal  far-field  pattemi  in  the  lateral  direction  can  be 
shaped  bv  varying  both  the  index  and  gain  profiles  (Fig  2<bi 
and  (cl)  Of  course  when  both  g(_v)and  n^vtart  symmetric, 
the  modal  patterns  are  symmetric  about  the  facet  normal, 
whether  n^v!  up  below  or  nsei  above  the  outside  values  a,  . 
a,|a,  •  a»  1  As  s  first  approximation  for  modeling  ttnpe  con¬ 
tact  lasers  we  take  a  •  a,  g.  •  g, .  a,  •  a,,  and  a.  *  a. 
Also  we  chooee  g.  snd  a.  ts  average  between  their  outside  and 
inside  valuer  t  e  .  g ,  •(/,-<»,  1,7  snd  a.  •  (m,  *  a,  )’2  Thus 
the  step  models  of  ff_v)and  n(_>-l  approximate  the  smooth  con¬ 
tinuous  distributions 

When  the  highgain  repon  is  centered  in  Is,  the  modal  pat 
terns  are  similar  to  the  H-G  functions  However,  as  repon  3  is 
shifted  padually  to  either  the  nght  or  left  of  center,  the  fat- 
field  patterns  become  asymmetnc  about  the  facet  normal 
The  mator  lobe  of  the  fundamental  mode  shifts  off  the  facet 
normal  For  example  when  A,  -  a,  -  10’*  and  when  repon  3 
®  shifted  off  center  by  I  urn.  the  ma;or  lobe  can  shift  almost 
10*  off  the  axis  normal,  and  the  higher  modes  have  loti  the 
characteristic  K-G  shapes  (I0| 

III  Comparison  with  Experimental  Patterns 
W>  studied  a  conventional  oxide  stripe  laser  7^77-9*4  whose 
patterns  llluttrate  the  asymmetry  typical  of  many  of  our  stnpe 
contact  lasers  Fig  I  is  the  cro«  section  of  the  starting  wafer 
336DH-CW.  which  was  pown  by  liquid-phase  epitaxy  The 
earner  concentrations  of  the  semiconducting  layers,  reading 
from  top  to  bottom,  are  10" /cm*  p  type.  3  *  I01'  p.  a  low 
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concentration  in  the  active  repon  from  a  trace  of  Si,  and  2  X 
IP1'  a  type  The  laser  had  an  oxide  stnpe  metallic  contact 
1 1 1]  about  10  sun  wide  and  cleaved  face  it  separated  100  urn 
One  facet  had  a  dielectric  reflecting  costing  According  to  the 
theory  of  current  spreading  in  a  two-layet  cap  (2),  the  effec¬ 
tive  width  of  the  recombtniuon  repon  u  10  urn  greater  than 
the  width  of  the  contact 

All  itudtes  used  a  low  duty  cycle  puiaed  current  (20  ns  nom 
inal  width.  5000ft)  to  avoid  significant  heating  Measure 
menu  were  made  with  *  box  car  integrator  with  10  ns  nominal 
gate  width  as  s  peak  detector  The  spontaneous  emission  was 
transmitted  through  a  tandem  pair  of  100  A  bandpass  filters 
with  centers  bOO  A  lets  than  the  coherent  band  Thu  short 
wave  emission  pves  an  account  of  the  junction  voltage  at  the 
output  facet  [2|  Further  details  of  the  technique  are  pven  in 
12) 

The  PI  tpuUn/auoft  index!  characteristic  was  nearh  linear 
with  soft  tum-on  No  marked  structure  was  present  m  It 
beyond  a  slight  wtvsnets  We  have  observed  all  the  effects  dis¬ 
cussed  in  lasers  with  very  lineat  characteristic  curvet  at  well  as 
those  with  sharp  kinks  tn  PI.  but  have  been  unable  to  correlate 
any  of  the  reported  properties  with  kinks  or  structure  The 
pattemt  here  presented  were  in  the  soft  tum-on  repon  close  to 
threshold 

The  model  psrsmeters  and  widths  of  repons  2.  3.  and  4  were 
admired  so  thar  the  theoretical  pattern  of  s  purr  mode  ap¬ 
proximated  both  the  monochromauc  far-field  pattern  and  the 
coherent  neat-field  The  luted  index  steps  are  consilient  with 
out  limned  knowledge  of  the  earner  concentration  beneath 
the  itnpe  contact  and  of  the  dependence  of  refracove  index 
on  earner  denim  |7]  With  the  gain  parameten  of  Table  I. 
the  modal  gain  g  •  -o  •  -  2  Real  (y)  u  about  1 00  cm ' 1  fot  each 
of  the  three  lowest  modes  At  threshold.  gm  (I  Tla,ra ,  s  con¬ 
dition  satisfied  by  T  •  0  6  and  nimf  ■  60  cm'1 

Fig  3  compsrei  the  expenmemal  near-field  patterns  with 
the  calculated  second -order  (S  •  2)  lateral  mode  The  shape  of 
the  spontaneous  emission  at  $0  mA  (threshold  *  40  mA)  u 
symmetrical  with  half  width  (full  width  at  half  power)  of  )4 
*rm  it  is  assumed  to  be  concentnc  with  hj  of  Fig  2  to  enable 
companion  with  the  calculations  The  coherent  pattern  at  $0 
mA.  P  *  4  mW.  (B  -ex  pen  men  tall  ts  a  single  lobe  with  faint 
structure,  half  width  of  8  srm.  and  offset  of  -2  urn  from  the 
spontaneous  The  calculated  pattern  (/>•  theory )  has  s  half 
width  of  3  8  son  and  displacement  of  -3.6  um  The  detail 
at  the  bottom  of  the  figure  locates  the  positrons  of  the  gain 
steps 

The  corresponding  beam  profiles  are  shown  in  Fig.  4  The 
solid  curve  u  measured  st  52  mA  and  wavelength  of  8192  A 
with  resolution  of  0  4  A  It  u  a  tingle  lobe  with  weak  minor 
structure,  half  width  of  7*.  and  peak  at  -6*  from  the  normal 
to  the  facet  The  theoretical  curve  (broken  line  I  is  similar  but 
with  slightly  less  structure,  and  has  the  same  half  width  and 
offset 

The  agreement  between  model  calculation  and  observation 
(tn  particular  of  the  beam  profiles,  where  the  instrumentation 
has  selected  s  pure  mode  rather  than  the  average  over  all  5 
modes  tn  the  entire  spectrum  svhlch  ts  recorded  for  the  near 
field)  gives  enough  credence  to  the  model  to  warrant  analyst 
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h$  J  W  Field  u)  Spontaneous  lasrsiuis  maanired  *1  ifcon  «r» 
taafth  lb)  Coherent  RIM  -  measured  and  calculated  (or  woM 
order  mod*  (c)  F  outsorts  of  steps  used  m  calculation  Canter  of 
term  reason  b4  of  Fra  2  assumed  to  be  on  center  of  spontaneous 
profile  Now  Out  Om  theoretical  and  aspertmeatsl  patterns  to  <b>  «r» 
both  shifted  to  Hint  dm c tron  Fitterot  of  leaet  7 ,'77. 4*4  10  mA 
•boo*  threshold.  R  •  4  ■« 

of  the  theoretical  pitttmt  of  different  5  modes  Fig  $  ihowt 
the  three  lowest,  with  detail  1 1 >  the  facet  illumination  and  (b) 
the  beam  profiles  One  itnkw|  feature  is  the  similarity  of  the 
shapes,  especially  m  the  beam  profiles  all  the  profiles  are 
dominated  by  a  single  lobe,  with  half  width  *arvtng  from  6  5* 
for  the  lowest  mode  to  4*  for  S  •  3  Abo  the  near  fields  of  S 
•  !  and  2  are  nearty  identical  and  only  by  5  *  3  is  pronounced 
structure  appearing.  Another  interesting  feature  of  both  sets 
of  patterns  n  the  direction  of  the  displacement,  which  changes 
agn  between  successive  modes,  in  the  figure  odd  modes  are 
displaced  to  the  nght.  even  to  the  left  The  calculated  modal 
dispersions  sre  AX ,,  *  *0.9  A  and  AX„  •  *2.5  A.  Fot  a  box 
mode  in  a  waveguide  of  width  14  urn  fthe  value  of  b,  in  Table 
I)  the  corresponding  values  would  be  -0.75  and  - 1 J  A  (I2| 

In  our  model  we  find  that  the  wavelength  and  the  near-  and 
far-field  patterns  of  each  mode  are  functions  of  the  dielectric 
profile  It  has  been  found  experimentally  that  the  gain  pro¬ 
file  above  threshold  ta  a  function  of  dnve  current  (2|.  (3| 
Hence,  the  parameters  of  Table  I  will  changr  with  current  and 
so  the  frequencies  and  patterns  will  changr  with  output  power 
Such  a  power  dependence  experimentally  found  ts  illustrated 
tn  Ftg  6  Detail  (a)  shows  i  high  resolution  (0  1 5  A)  spectrum 
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Fig  4  hears  profiles,  conditions  of  Fg  )  Both  paiiemi  hiv*  a  ungk 
lobe.  shift  above  -  J* 


*os*ov  04  tact'  <»ai  nMVs.  aaasl  Mt 
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F«  S  Noar  an*  fsrfseld  calculations  three  lowest  modes  Note  dial 
beam  ptoflWi  has*  only  a  waste  lob*  la)  Facet  illumination  (bt 
Beam  proflk 


efcMVO  ill  OH>tc*«s»  im* 

ta)  (b> 


Fig  i  Section  of  a*ectfiHi  and  beam  profile  ai  different  currants 
Nose  the  needs  ehlfi  with  current  of  podtson  of  peak  Detail  la)  her 
one  norm aiuj non  detail  (b)  another  Rower  under  SO  mA  trace  • 
0  2  mA  Lamt  777-S.( 

of  the  beam  collected  by  a  large  aperture  aspheric  lent  for 
total  power  ran  png  from  3  to  6  mW  (the  power  undet  the 
doited  line  for  50  mA  is  about  0.2  mW>  Nonce  the  shift  with 
current  of  the  wavelength  of  the  mode  with  no  evidence  of 
excitation  of  other  lateral  modes  until  65  mA  A  cunous  fee- 
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tuit  no:  predicted  b>  the  model  m  iu  present  soptusucauon  u 
the  wide  line,  shout  0.5  A. 

Mb i  gjvet  the  corresponding  Uteril  profiles  The  u 
phene  lens  has  been  replaced  with  a  cylindrical  one.  the  slat 
widened  to  a  resolution  of  0  4  K.  and  the  spectrometer  tuned 
to  maximum  ugna:  at  each  current  Nonce  the  nearly  perfect 
repetition  of  the  shape  I  the  curve  for  65  mA  duplicates  that 
for  60  mA  except  for  an  extra  shift  of  2*  and  slifht  change  ui 
power )  The  one  obvious  change  with  current  is  in  the  poo 
non  of  the  peak,  which  seems  to  shift  monotomcallv  from  *  2* 
at  *8  mA  to  -9*  at  65 

Experimental  features  of  the  many  stripe  lasers  'with  asym¬ 
metrical  beams  that  we  have  studied  are  the  shift  of  wave¬ 
length  and  field  profiles  with  current  and  the  broad  spectral 
lines  Beam  and  spectral  shifts  with  current  can  be  ued  to 
the  current  dependence  of  the  dielectric  profile  Although 
it  is  conceivable  that  the  Unewidth  is  i  measure  of  the  change 
in  current  during  the  10  iu  sampling  gate,  the  required  shift 
during  the  sampling  period  seems  excessive  to  account  for  the 
broad  lines  Moreover  the  same  Unewidth  is  observed  unde; 
CM'  operaoor. 

rv  DtJccssiON 

In  out  step  model,  important  asymmetry  m  the  radiation 
occurs  only  when  there  u  ipprecubie  detocusng  by  the  index 
profile  Modes  with  n t  >  n,.  n,  produce  very  little  asym¬ 
metry  Only  with  s  large  neptm  step  in  index  plut  gain 
guiding  did  we  find  asymmetries  of  the  sire  observed.  The 
value  oaed  in  the  calculation  An  ■  -  4  X  10*’  teems  to  be 
compatible  with  our  present  undemanding  of  the  perturba¬ 
tion  tv  the  index  accompanying  iniecnon  of  enough  free  ear¬ 
ners  to  give  the  required  locau  gain  of  300  cm*1 

Because  of  the  near  cancellation  of  gam  guiding  by  index 
anuguiding.  both  of  which  change  with  mtecnon  level,  the 
propagation  is  quite  lensitive  to  change  m  gam.  which  makes 
the  propapbon  constant  i  function  of  current  tor  power) 
This  rubberv  nature  of  the  cavity  mav  be  the  cause  of  the  large 
changes  with  current  of  ’die  modal  wavelengths  and  the  non¬ 
monotonic  current  dependence  of  modal  power  ihowr.  in 
Fig  6 

It  is  difficult  to  deads  the  efficacy  of  replacing  the  actual 
dieiectnc  and  gain  distribution  in  the  later  with  the  step  ap¬ 
proximation  When  the  step  3  is  centered  in  the  guide,  re¬ 
moving  the  Mymmetry.  the  modal  patterns  become  nnular 
to  H-G  modes,  which  seems  to  ipeak  for  the  utility  of  the 
model  A  poaubte  cause  of  the  asymmetries  introduced  m  Fig 
2  is  the  imperfections  m  the  stripe  contact  Poaible  sources  of 
contact  asvmmetrv  art  shadowing  r.  opening  the  window 
through  the  oxide  layer  or  in  dcpaauon  of  the  metal  While 
we  have  modeled  the  datornon  as  a  Axed  displacement  of  the 
extrema  of  both  gain  and  index,  it  seems  probable  that  da- 
placemen:  of  one  relative  to  the  other  would  wrve  M  well 
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The  latter  has  the  attractive  feature  that  the  distortion  need 
not  affect  the  current  distribution,  agreeing  with  our  observa¬ 
tion  of  a  quite  symmetric  spontaneous  pattern  on  the  facet 
accompanying  strongly  distorted  coherent  patterns. 

In  conclusion  we  have  shown  with  a  specific  example  that  an 
asymmetric  siep  waveguide  can  produce  the  type  of  distorted 
patterns  observed  in  many  itnpe  lasers,  that  for  it  the  depen¬ 
dence  of  partem  shape  on  mode  number  ts  so  weak  that  the 
modal  content  cannot  be  readily  estimated  from  the  beam  pro¬ 
file.  and  that  when  index  defocuang  u  comparable  with  gain 
guiding  the  propagation  constant  becomes  very  sanative  to 
changes  in  gain. 


ACKNOWLEDGMENT 

Wafer  536DH-C  W  was  supplied  by  M  Ellen  berg  We  also 
acknowledge  the  excellent  technical  assistance  of  L  EJbaum. 
M  G  Harvey.  A.  kart.  H  V  Kowgtt.  and  D  P  Manneili 

References 

I  I  T  K  Zacnoi  an <3  J  E  Rippe:  "Resonant  model  of  GaAi  iuiu 
uoo  lasers."  IEEE  J  ffttmmm  EUem>K.  vol  QE-!.  pp  29-37. 
Jen  1 969 

[2!  H  S  Sommers  it  end  D  O  North.  "Ex  penmen  let  end  theo- 
reucel  nudv  of  the  mens!  venation  of  junction  voltage  end  cur¬ 
rent  dietnbunon  in  narrow  unpe  injection  leaers  *  J  Appi  Art 
voi  a*  pp  44*0- aod’  Oct  i*rr 

( 3 1  F  A.  KirSby.  A  R  Goodman.  G  H  I  Thom  peon,  and  P  R 
Setwai.  "Obeemoom  of  relHocuang  in  imps  geometry  ttnu 
conductor  laeert  and  the  development  of  a  comprehennve  model 
of  their  operation  IE  El  J  frimmm  Drrm .  vol  QE-I }.  pp 
70S-' 19.  A  up  197'.  K.  Kobavaahi.  R_  Lang  H.  Yoottu.  L  ialima 
and  1  Haveslu.  “Homouiai  mode  deformation  and  anomaloui 
taamg  properuei  of  itnpe  goometn  injection  laaoivta  pan  mental. ~ 
Jspm  J  Ap pi  Ah vi .  vol  |g.  pp  207-20*.  Jen.  1977 
(4]  N  Qunont  “Nonbneemv  ui  power -output-current  chars  cams 
tics  o!  irnpe-peometrs  injection  laaen  -  J  Appi  Art  vol  4| 
pp  3237-324}  a  up  197'.  introduces  an  asymmetry  »  gnm  pro¬ 
file  ai  perturbation  to  the  H-G  treatment,  and  etammes  the 
•fleer  or  the  fundamental  mode  W*  are  conamed  with  the  dis¬ 
torted  *<r*-om rt  mode*  which  are  escsted  m  our  narrow  ttnpe 
laaen 

1 3  j  D  D  CooS  end  F  R  Netfi.  Gain-induced  ftudmi  and  aaugmauc 
output  beam*  of  G*Aj  laaen  ~  J  A  pot  Art .  vol  a*,  pp  14*0- 

1*72.  Apt  1 97! 

I*i  I  W  Hakkx  “Smped  G*Ai  linn  Mode  nu  and  tfTicatnrv .*  J 
Appi  Art  vol  «*.  pp  2723- 2“ 30  June  197! 

(7)  T  L  Ho*.  "Weveguadma  m  a  rtnpa  |aumetr>  tunctaon  ln«i.~ 
IEEE  J  Onenrm  Cermm.  vol  (JE-lJ  pp  UJ-**i  Aug 
1977 

fit  E.A.J  Maianh  DlMectfSc rectangular wpvugmde mm dnocooMl 

coupler  fa  enegreied  opnea.-  M  Jen  Teed  y.  vol  aa.  pp 

2071-2102.  Sept  19*9 

[9|  J  1  better  and  H  KraarnL  “Dtngn  curves  tor  doubts  beserw- 
Juwcooa  Ml  drodaa."  MCA  Mn  eel  31.  pp  MJ-SJl  Don 

1977 

{10)  L  I  Setoff.  Ootm me.  MrcAemci  Not  York  McGnwdfBL 

1935.  p  *0 

[111  l  Ledenv  and  H.  Kremel.  The  influence  of  device  fabnaooe 

dtodUr^appr  Avr  tJn',,wM,^J^T0^7i^,CS^,i j^pTa” 
[12)  H  S  Sommers  Jr  and  D  0  North  “Tba  power  pinne  of  m- 
jecooe  leans  The  theory  and  aepawmmt  pp  a  nonlmetr  model 
o <  Man*."  JpdM-Jwie  Clrcmm  .eel  IS.  pp  *73-499.  197* 


II1EI  JOURNAL  Ol  QUANTUM  ELECTRONICS.  VOL  QI-14.  NO  *>.  JULY  !•*»• 


*0** 


A  Rigorous  Boundary  Value  Solution  for  the  Lateral 
Modes  of  Stripe  Geometry  Injection  Lasers 

JEROME  k  BITLER  MU4MK  ieei  and  JOSEPH  B  DELANEY 


■INWfl  -  A  MM  nfkwwslKil  noM  oootol  las  awitysac  lataral 
lotos  a I  uiyt  nwim  Iwn  »  pimM  TW  isilto  iinfi  Uan  kk 
■  o4*lsl  a  •  tks**-l*y*s  an  n»rii  la  wfcick  iIh  toatertstc  ctutnn  of 
tto  mm  lira  mm  owb  siosif  to*  totesal  dssvcuow.  to*  intolWi 
immti  ol  tto  Nm«Ma|  pan  taym  ■  HMMd  to  to  MHM 
totfMtoH  Tto  stoaooa  torOotpin  tflorto  •  nfutow  ■  lift  mg  ol 
toe  ftoto  ol  to*  Kim  lorn  otto  too**  of  ito  swssiwitiwg  fto" 
Inm  To  ttoMtni*  il»  motor  dw  aotot  of  s  nwftdi  vitk  pan 
Mu  ato»rss».  nnaim  atuwf  dw  Utasal  dimm  at  muigimt 
Tto  fwlto  on  »niw*  a  •  tawai  ctaatoaulna  of  Waaia-fawaa 
(M-f.i  (ucuowt.  toniofw*  ttoto  km  tom  Ssscrttod  witfc  •  waft* 
H-C  foactut  fwfanaul  awad*  *lal  ttoot  an  ti  •Mrowvii 
a  cakatoiod  tad  ratoito  to  dw  pm  dwtnbauoa  (Psotsums  nnautot 
ol  dw  kanl  ft*M  and  of  dw  faadnaul  aiod*  auag  •  ua#« 
H-C  faacnua  aot  npraak  ailcM  at  Ito  kn— dsswi  com  yMta 
tom  am  w  aack  a  (0  pawl  toffmai  (toai  maiu  catcatoMd 
from  latm  naatoa»l><w  of  H-C  faacltoat  I  la  addiooa.  dw  pul 
$oto  fWdi  am  to— M—d  at  tkmdtoM  fat  nssom  aart  swsds 
fooooovnm 


iNTRODl'CTION 

THF  modal  characteristics  of  ttrtpe  geometry  infection 
laseri  are  preiented  here  The  theoretical  mi-del  de¬ 
veloped  allows  for  a  rigor  out  treatment  of  both  lateral  (along 
the  junction  plane  i  and  vertical  modes  of  the  oxide  ttnpe  in 
pection  latert  Earh  work  focused  on  a  mode!  which  treated 
the  active  region  at  a  light  guiding  ttructure  with  a  imocth 
cor.tinuout  dielectric  variation  along  the  transverse  direction! 
at  opposed  to  abrupt  rndra  tlept  In  particular  the  aatumed 
dielectric  variation  along  the  lateral  (a  direction)  and  the 
vertical  (i  direction i  wai  parabolic  of  the  form  six.  vl* 
toiMX1  *  by*  I  1 1 )  The  accuracy  of  thn  detcnption  it  not 
unreasonable  unce  it  wai  applied  to  modeling  homoiunction 
laaert  However  the  model  n  not  applicable  to  heteroiunction 
devicet  where  the  grown  active  layer,  generally  very  thin  for 
contemporary  Ok'  laaert.  it  sandwiched  between  two  passive 
laven  with  smaller  dielectric  constants  |2)  A  more  accurate 
detcnption  of  the  active  layer  it  one  with  a  refractive  mdea 
which  is  constant  along  the  vertical  direction  hut  vanes 
parabolically  along  the  lateral  dimension  [T] .  |4|  This  di¬ 
electric  detcnption  it  reasonable  since  earner  m version  u 
uniform  along  the  vertical  dimenuon.  whereat  the  inversion 
a  not  uniform  along  the  lateral  dimension  due  to  current 
spreading  from  the  stripe  contact 
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Standard  techniques  for  solving  boundary  value  problems 
do  not  readily  apply  to  the  geometry  of  the  regions  ui  the 
vicinity  of  the  active  layer  Consequently,  prenout  analyses 
of  the  lateral  mode  ttructure  of  contemporary  Ok  stripe 
geometry  devices  have  been  approximate  In  particulai  the 
dielectric  constant  of  tht  active  layer  vanes  along  the  lateral 
dimenuon  while  the  dielec'ric  constant  of  the  bounding  p-  and 
n-AJCaAs  layers  it  homogeneous  For  parabolic  index  vans 
lions  in  the  active  layer,  the  solution  of  the  wave  equation 
yields  Hermite -Gaussian  (H-C I  functions  On  the  other  hand, 
solutions  of  the  wave  equation  in  the  taven  bounding  the 
active  region!  are  not  H-C  functions  Appropriate  field  quan 
titles  cannot  be  matched  at  the  boundaries  The  method  used 
here  incorporates  a  technique  developed  previously  |5)  The 
solution  of  the  w»vr  equation  tn  the  active  layer  is  written  as 
a  linear  combination  of  H-C  functions  whereas,  in  the  bound 
mg  patuvr  regions  the  solution  of  the  wave  equation  is  writ¬ 
ten  as  a  linear  combination  of  plane  waves  Using  boundary 
conditions  at  the  heterotunctiont.  an  eigenequation  it  devel¬ 
oped  whose  solution  gives  the  appropriate  modal  propagation 
constant!  The  normal  modes  of  rhe  laser  structure  are  linear 
combinations  of  H-C  functions  with  weighting  coefficients 
dependent  upon  the  specifics  of  tise  geometry  The  theory  is 
developed  to  the  extent  that  complex  dielectric  constants  can 
be  used  in  the  model  This  meant  that  itructures  which  have 
either  lateral  index  ot  gam  guiding  properties  can  be  analyzed 

An  understanding  of  lateral  mode  stability  of  stnpe  geom¬ 
etry  lasers  has  led  to  interest  in  accuratr  modeling  of  the  lat¬ 
eral  index  profile  and  of  lateral  mode  characteristics  Lateral 
mode  perturbation  due  to  nonparabolicitv  of  the  gam  profile 
was  also  treated  by  expanding  the  fundamental  mode  in  a 
senes  of  H-C  functions  |b)  A  complete  model  fot  the  lateral 
mode  characteristics  would  then  include  these  perturbed 
solutions  due  to  aberrations  of  the  parabolic  gam  profile 
property  matched  at  die  boundaries  as  we  have  discussed 
However,  in  Uus  paper  we  address  ourselves  solely  to  die 
effects  on  lateral  modes  of  matching  solutions  at  die  bound 
anet  and  'Uustrate  these  effects  using  a  parabolic  dielectnc 
model 

Theory 

Due  to  the  nature  of  the  geometrical  structure  of  stnpe 
geometry  injection  lasers.  It  ts  difficult  to  formulate  a  simple 
eigenvalue  problem  accurately  characterizing  the  modes  As 
mentioned  previously .  the  eigenmodes  of  stnpe  geometry 
taarn  can  be  described  as  linear  combinations  of  H-C  func 
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uam  obtained  from  t  solution  of  Maxwell's  equations  The 
field  equations  in  reduced  to  in  approximate  form  wtuch 
Oescnbei  the  waveguide  model  having  (tec trie  field!  poianied 
cither  aiong  x  ot  y  as  shown  tr.  Fit  1  The  fieidi  with  x  and  y 
putaruatioiu  ut  defined  by  the  wti  of  model  T}m  and 
respectively  [T]  Since  moit  iniection  iisen  have 
fields  poLarued  predominantly  alon|  x  (S|  we  Umi!  our 
discussion  to  the  T  *  model  Aa>umin|  mode  propagation  ha> 
the  form  expiito/  -  vii  where  y  •  -C  Z  *  C  the  field  equa 
tioni  become  (the  superscript  1  u  dropped  1 


£,  ■  ♦tx.y) 

(1) 

£,  •  0 

(2) 

_  1  J* 
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„  ._L_  il* 
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H  --  1  Iri. 
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t,n,  to 

where  k«  is  the  free 
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71.  and 

atx.  vi  is  the  relative  dielectric  constant  The  transverse  wave 

function  +  lauafiet 

7 f (7) 

In  the  ncmiry  of  the  actree  region .  the  waveguide  geometry 
hat  a  form  illustrated  by  Fig  2.  In  the  actm  layer,  lateral 
field  confinement  ts  due  to  the  following 
11  the  nonuruform  distribution  of  the  tndea  of  refraction 
and 

2)  the  nonuruform  distribution  of  the  fam  losses  [9] .  ||0] 
A  re  aeon  that  maxes  determination  of  the  mam  tourer  of 
lateral  confinement  difficult  ts  the  fact  that  the  optical  (am 
lots  coefflcwnt  and  the  index  of  refraction  are  inextncably 
related  On  the  other  hand,  the  source  of  optical  confinement 
along  the  vertical  v  direction  can  be  eauh  determined  The 
large  refractive  index  diacor.imuities  at  the  heterotunctions 
which  bound  the  active  layer  are  aimost  solely  responsible  for 
both  optical  and  earner  confinement  to  the  actm  layer 
To  formulate  the  Laser  modes  the  waveguide  wUl  he  charac 
tamed  as  follows  The  taro  horuontaJ  petal ve  liven  A  and  C 
have  complex  dielectric  constant!  which  are  petition  aide 


pendant  The  refractive  index  of  them  layers  is  determined  by 
the  lasing  wavelength  and  aluminum  concentration,  the  bulk 
absorption  coefficient  is  positive  The  active  layer  B  ts  assumed 
to  hive  a  complex  disiectnc  constant  which  a  dependent 
only  on  the  lateral  1  direction 

In  the  following  mathemaucal  formulation  we  assume  that 
the  dielectric  constant  u  rvmmetncal  about  x  ■  0  In  the 
active  layer  the  wave  equation  has  the  form 

7f*»v  >’)♦,•  0  (8) 

and  six  1  u  the  complex  dielectric  constant  of  the  active  layer 
which  has  only  lateral  variations  In  the  active  layer,  the 
solution  to  (8)  can  be  obtained  by  separation  of  variables .  we 
write  a  typical  solution  of  the  form 

♦*U.v)»  (Mxldafy)  (9) 

where  the  vertical  function  •(.  u 

d*(y)  •  cot  qy  (10) 

The  lateral  field  function  satisfies 

♦  (kiaul-g1  vy^ba-O  (II) 

ax 

As  u  weO  known,  there  exist  rwo  sets  of  vertical  modes  b»(yl 
where  ( 10*  n  the  solution  of  even  modes  about  y  •  0  The 
modes  can  have  discrete  g  values  representing  trapped  modes 
or  it  can  take  on  a  range  of  continuous  g  values,  radiation 
modes  along  the  y  direction  The  number  of  trapped  vertical 
modes  a  dependent  upon  the  cavity  width  d  and  the  index 
steps  between  the  active  region  and  neighboring  regions  A  and 
C  fg) .  Contemporary  O  injection  lawn  have  narrow  active 
layers,  d  —  0.1 -0J  ion.  and  moderate  index  steps  Them 
structures  generally  support  only  the  fundamental  trapped 
mode  Consequently  we  limit  our  discussion  here  to  only 
the  fundamental  vertical  mode 

We  now  conndet  the  solutions  of  ( 1 1 1  which  give  the  lateral 
mode  character  It  should  be  noted  that  (11)  hat  1  one- 
dimensional  form  arrulai  to  typical  slab  waveguides  For 
ciampfe .  Use  quantity  y*  -  q1  can  be  interpreted  as  1  one 
dimen wonal  eigenvalue 

The  solution  for  b»U)  can  be  effected  once  the  value  of 
■  (xl  is  specified  (n  terms  of  the  bulk  parameters  <Kxt  the 
local  refractive  index  and  eul/gui  the  local  absorption/ gam 
coefficient  the  dielectric  constant  can  be  approximated  as 

atx)  *  n*(x)  •  i«u  hnxVka 
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In  the  cjje  of  negative  abiorpdon  or  gam  vat  write  umpl> 
at*)*  f^x)  In  Fie  3  we  illustrate  two  index  absorption 
gain  profiles  In  Fig  .Mai.  the  index,  fair,  has  a  paiabolic  pro 
file  while  in  (bl  the  optical  parameters  arc  shown  to  have  step 
changes  Before  we  discuis  the  details  of  a  specific  structure 
we  note  that  the  lateral  modes  obtained  from  the  parabolu 
index  mode!  are  only  discrete  trapped  modes  whereas  the 
step  index  model  has  a  set  of  tiapped  as  well  as  radiation 
modes  In  our  discussion  onl)  the  trapped  modes  will  be 
used  and  for  the  sake  of  convenience  we  ..onsider  the  para 
bolic  index  ptoftit  Therefore  the  dielectric  constant  in  the 
active  lavet  is 

ala) •  a,  -  * ;«’ t3  (13) 

where  a„  is  the  complex  dielectric  at  »  *  0  From  1 12).  one 

l««» 

a,  ■  "!(0t  ♦  iflOWOi  l*0  I04l 

l  sing  the  parameters  "eVRtOl.go  •g'O  ft"  if.  and  x0  ** 
have 

»i  3  i2"0ft" 4 1 "a*/  **  10* yu»*o»! 

where  x0  is  in  urn,  *„  has  dimensions  urn'1  .g  Sf  are  in  cm '' 
The  wave  equation  becomes 

</3w» 

~r  *  !*««•  *  t*  •  «8*i*,|  v>  •  o  d*i 

ax 

The  solution  of  <I4|  yields  the  standard  H-G  wavefunettoni 
However,  if  we  take  a  single  H-G  function,  then  the  boundary 
conditions  at  the  hcteroiunction  interlaces  cannot  be  met  It 
is  therefore  necessary  to  define  a  mode  as  one  being  a  linear 
combination  of  H-G  functions  The  solution  has  the  form 
(the  b  subscript  on  wli  l  and  Mil  is  temporarily  dropped  for 
notations!  convenience) 


In  ( Ibbl  H,  is  the  Hcrmite  polynomial,  and  q,  satisfies 

*«*<•  *7S  *  q}  *ks,M2l*  1)  (17) 

The  lateral  modes  can  be  divided  into  two  sets  depending  upon 
whether  the  fields  are  even  or  odd  about  x  •  0  For  even  fields 
the  index  /  ■  0.2  *  whereas  the  odd  fields  have  /  *  1 .  3. 
5.  Ike  limit  out  discussion  here  only  to  the  even  fields 
since  primary  interest  will  be  focused  on  the  fundamental 
waveguide  mode 

A  transverse  mode  is  of  course  defined  by  a  given  propa¬ 
gation  constant  y  Therefore  ♦„(* ,  y )  in  1 1 5 1  gives  the  trans 
verse  modal  field  distribution  with  the  total  field  solution 
being 

♦  »(*.. s  te"’*  •  e'v  21  A,q,(y)  w/U)  <  1 8 1 

#•0 

We  see  that  the  right-hand  side  of  (Pi  must  be  independent 
of  the  index  l  fcven  though  q,  the  vertical  "wavenumber” 
varies  with  /.  the  vertical  y  dependence  in  (15)  defines  a 
specific  waveguide  mode  which  in  out  case  is  the  fundamen 
ta!  vertical  mode  This  particular  aspect  of  the  problem  will 
he  better  understood  as  we  develop  the  total  transverse  mode 
eigenvalue  solutions 

For  the  sake  of  simplicity .  we  considei  a  symmetric  wave¬ 
guide  structure  along  the  vertical  so  that  +4(x  y  I  ■  4-,.(x  -  vl 
•  double  heterostructure  GaAs-AIGaAs  lasers  are  generally 
fabricated  with  equal  amounts  of  aluminum  in  the  p  and 
n-walls  sandwiching  the  active  layer)  Outside  the  active  layer, 
the  w  avefunctson 


F,(.t  vt»  f 


#(X>  cos  x1  exp 


(tv)(x3-*; 


which  is  symmetric  about  x  •  0 
To  determine  B(\  \  and  the  propagation  constant  y.it  is  neces 
sarx  to  applv  the  appropriate  boundary  conditions  at  the 
hetero)unctions  Therefor*  at  .1  •  d'2 

21  |  B[\)cm\xd\  (20) 

From  the  inverse  Fourier  transform  £h\l  is 

L'Mr(!)*#(X)  <211 

where  we  have  defined  the  transform  pair 
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The  matrix  elements 
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or3  •  *;  •  d;ix>*x  (2*) 

The  propaption  coniurt  y  is  now  determined  from  (2bi 
The  nontrivial  iciutton  a  pv»-  u  root!  of  the  equation 

(ktt£-n>»o  i:s> 

where  £  u  the  unit  nutm  Note  that  tin  matrices  E  and  n 
have  infinite  otde-  Consequently ,  it  u  impossible  to  find  the 
eigenvalue*  '  s%  itandarc  techniques  One  approach  a  to 
u*ume  i  I  x  1  maim  and  determine  y  front  the  retultinf 
transcendents  equation  Next  expend  £  and  £2  to  2  x  I.  etc 
The  resulting  tequence  of  y  value*  thus  forms  t  Cauchy  se¬ 
quence  If  the  ttquence  spproache*  i  'muting  value  at  the 
mat  rut  orden  increase  then  the  sequence  converges  In  Table 
I  we  show  the  fundamental  mode  y  value*  for  different  matm 
order*  Note  alto  the  expansion  coefficient!  A,.  A).  A, 
can  be  calculated  after  the  appropriate  eigenvalues  an  found 

For  the  tace  of  c ample tenet! .  at  thn  pom i  we  reconsider  the 
vanoui  modes  involved  in  this  formulation  Different  lateral 
mode*  are  found  from  solutions  of  (28 1  For  example  for  a 
I  x  [  matrix  there  enll  Sc  only  one  eigenvalue  the  one  *p- 
peoxsmaunf  the  fundamental  mode,  while  for  t  2  x  2  matrix, 
there  exist  two  eqten  values  those  of  the  fundamental  and 
first-order  mode*  As  the  matrix  orden  increase  the  com¬ 
plete  set  of  lateral  modes  is  calculated  The  vertical  mode* 
can  be  found  from  (281  by  judicious  choice*  of  q  To  under- 
itind  thn  consider  fint  the  simple  1  *  1  equation  with  H-G 
function! 


*»(^7r)«P  (-~)  (»> 

so  tha: 

“•(••t)  f 

<f  •  *u’|^»7j  • 

exp^-^-jdj  (301 

In  the  Umitui|  case  a  —  0  mo  curvature  I. 
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The  resultuii  eiftncquation  D««  *  I  u  identical  to  that  of 
the  three -la  vet  slab  wgveguide  where  p,  •(-y1  -  k;)in  <sA 
pi  ♦  qj  “*«(**  *  ».l  Consequently,  vertical  modei  found 
from  (31)  can  be  characterued  bv  different  q,  vaiuei  found 
from  pe  ■  de  '-an  idt>d.‘2)  The  q,(l  ■  2.  *.  i  value*  uaed  in 
(28 1  art  very  dote  to  qe  wher.  a  u  small  t  typical  itnpe  laatr 
structure*!  From  ( lTf  ere  fmd  qf  ■  q\  •  lie 

Numerical  Solutioni 

Several  numerical  tolutiom  pertaining  to  contemporarv 
itnpe  feometn  laaeri  will  be  presented  now  In  particular  ere 
will  study  the  convergence  of  the  eigenvalue*  pertain m|  to  the 
fundamental  eraveguide  mode  and  fundamental  mode  spot 
axes  a*  a  function  of  device  geometry 

Aa  mentioned  prmouah  the  solution  of  the  eigenvalue 
equation  (28)  cannot  be  obtained  using  matrices  of  infinite 
order  Therefore,  to  determine  the  trend  of  the  fundamental 
mode  field  solution*  we  first  assume  that  the  matrix  orden  of 
£  and  O  art  I  x  1  and  increase  their  orden  up  to  4  x  a  and 
observe  convergence  of  fields  md  their  eigenvalues  The 
specific  parameter!  used  m  our  model  defining  the  parabolic 
dielectric  constant  are  given  ui  Table  I  In  the  table  ,  we  show 
the  values  of  the  normalised  propagation  constant  for  the 
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various  orders  and  the  expinuon  coefficient!  A,  where  it  it 
arbitrarily  assumed  that  A,  ■  !  ♦  10  for  ili  cases  Sole  that 
for  the  four  component  case  the  magnitude  0f  the  coefficient! 
decreases  rapidly  which  me  am  that  the  ttructute  ir  region  B 
hat  an  a  dependence  of  the  form  of  a  Gaussian  function 
These  expannor  coefficient!  are  ilrongly  dependent  upor. 
geometrical  parameters  such  ai  the  cavity  width  J  and  dieie. 
trie  ttepi  The  rate  of  convergence  of  the  A ,'s  slow*  at  d  de 
crtairi  Thu  meant  that  for  thin  cavities  high  percentages  of 
high  order  M-G  functioni  are  preient  in  the  fundamental  cjv 
lt>  mode  It  thouid  he  noted  that  the  value  of  t  changes 
with  matrix  order  however  the  difference  between  luccei 
uvr  values  decreases  with  order  In  general  it  u  found  that 
the  single  -component  cate  n  reasonablv  accurate  for  it  rue 
turei  which  have  iargr  J  valuei  or  wher.  the  optical  field 
energy  t$  confined  to  the  active  report  b 

The  lateral  fieid  dittnbution  of  the  fundamental  »a»eguide 
mode  ii  shown  in  Fig  4  Curve  C ,  repreienti  the  lateral  field 
profile  when  the  matrix  order  it  1  x  I  tone  component  t  while 
C«  pertami  to  the  matrix  of  order  axa  tfour  component  i  It 
a  important  to  note  that  the  curve*  Ci  and  C,  are  almost 
identical  w-fuch  impliet  that  the  eigemaluei  and  fietdi  are 
converging  rapidlv  with  respect  to  matrix  order  It  u  seer 
that  the  effect  of  the  high -orde-  termi  of  the  field  expansion 
is  to  mcreaae  the  lateral  distribution  of  the  field  or  to  increaw 
the  field  spot  sue  Ptivscallv  the  spot  sue  increase  results 
from  the  optical  field  m  repon  fi  leaking  energy  into  the 
surrounding  lasers  A  and  C  Thu  lateral  spreading  effect  is 
due  to  the  fact  that  there  is  no  mechanism  to  lateralis  confine 
the  waveguide  modes  in  regions  A  and  C  lateral  confinement 
is  due  solely  to  the  graded  dielectric  constant  in  the  active 
layer  B 

the  next  turn  our  attention  to  calculation!  of  the  lateral  spot 
sues  at  halfpowtr  as  a  function  of  the  gam  distribution  param 
eter  2x,  In  these  calculations  we  define  x»  as  the  pom: 
where  Lateral  gam  is  aero,  i  e  .  &g  •  g,  The  value  of  g0  will  hr 
appropriately  adiusted  ao  that  the  modal  pm  C  •  2Real 
!>»•  SO  cm"'  which  accounts  fo*  radiation  end  losaes  This 
value  approximates  the  end  loxaes  of  a  later  of  length  /  •  240 
utt.  In  the  various  curves  obtained  from  numerical  calcuta 
lions,  we  assume  that  the  mdex  of  refraction  at  a  given  pom: 


ft  i 


along  the  latcial  direction  is  a  tuncuon  of  the  gain  or  earner 
inversion  at  that  point  When  the  value  bn  n  negative  the 
mdex  of  retraction  piovidei  an  amiguidtng  effect,  while  if  bn 
is  positive  the  refractive  index  provides  a  guiding  effect  In 
these  calculation!  we  assume  the  relation 


A  *  6n/6g  cm  (32) 

where  the  parameter  A  represent!  the  effect  of  the  put  on 
the  mdex  Wher.  A  ■  0  the  index  hai  no  curvature  Assume 
for  example  bf  *  100  cm"1  and  A  •  10"'  cm  then  bn  •  0.001 
In  Figi  5  and  6  we  show  the  total  width  of  the  fundamental 
mode  at  halfpower  as  a  function  of  the  gain  width  2x,  for 
index  steps  n„  n,  •  0  !  (Fig  5l  and  0.2  (Fig  6).  in  (a|.  the 
total  cavity  thickness  is  02  tarn  while  in  (bl  the  width  is  0  I 
gam 

l  sing  the  single  H-G  function  to  describe  the  fundamental 
mode  the  spot  sue  A  n  pver  by  |3) 


„  2  (In  2),n 

/j  ■  - . - 

I" 1,4 

where  T  is  the  vertical  field  confinement  factor  and 

e.  •  Reallaf  ("*** *  -  Rea!  (2A  ♦  i  I0"4/*,)*'3 
4.x, 

The  peak  pm  g,  »  gr*er  by 

.  P  4  II  «  FkAxSn.GlO-*)1”)3  4  1 _ T 

,e  ’  R4.xin»  10"*  *  f  °4 
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(34 1 


(35i 


The  term  (I  I'm,  T  in  ( 3 5 »  represents  lots  in  the  p-  and 
n-repons  adiacent  to  the  active  lavet  when  o,  •  a,  |I3| 

As  an  example  lake  d  *  0.2  um.  2x,  •  10  urr.  X  •  0  urn 
G  *  SO  cm"'  a.  •  a,.  *  20  cm"' .  it.  ■  3  6.  n,  ■  n,  •  3  4,  and 
A  *  0  The  confinement  factor  T  •  0  6  so  that  (33 1  pves  B  • 
"  Jb  um  From  our  more  exact  model  we  find  a  spot  sue  of 
F  02  um  an  F  percent  difference  In  genera!  we  find  that  the 
computed  spot  sue  uung  our  model  differs  from  the  approxi 
mate  B  by  n  much  as  30  percent  »e  show  in  Fig  Mdashed 
curvesi  the  approximate  beam  sue  B  found  from(33l  The 
relative  comparison!  between  the  approximate  curvet  and  the 
numerical  ones  are  smtsia-  fot  other  d  values  and  index  steps 
The  active  peak  gam  Ic  is  shown  in  F-P  T  and  F  We  note 
igiit!  that  g,  is  chose*  such  that  the  modal  pm  is  50  cm"1 
to  account  for  end  losaes  The  two  sell  of  curves  ate  fot 
d»  0  i  and  0.2  aim  the  required  gain  for  each  sei  starts  to 
increase  for  2x.  valuei  below  -  10  wm  Finally  sve  show  m 
Fig  ®  the  pep  pm  g ,  at  a  function  of  cavity  thickness  D  for 
a  given  pm  region  thickness  2x. 


Cosset  r»on 

A  new  mathematical  model  describing  the  waveguide  modes 
of  the  stripe  protnetry  in  lection  case:  hat  been  presented 
Previous  models  usrd  to  describe  these  lasers  have  not  giver, 
sufficient  detail  to  the  interaction  between  the  vertical  and 
latrral  fields  and  to  the  nature  of  the  boundary  between  the 
acme  lave-  and  the  paasrvr  p  and  n-AIGaAs  walls  The  mode! 
prewntrd  here  was  applied  to  structures  with  dielectric  step 
variation!  m  the  vertical  direction  and  with  petabolu  vana- 
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non*  of  th*  dKMcinc  contrat  iloog  th*  lettrai  dir*ctton 
lateral  vtnatKm*  war*  umited  to  the  K!t»*  1 *var  Although  w* 
it*  tL  1 1 tad  th*  mod*  ttmcrurt  of  th*  Imtti  parabolic  *ary  mg 
tehctnc  th*  modal  can  b*  aa ad  to  anafyw  arbitrary  domett* 

To  and* nt and  th*  nature  of  th*  ft*M  aotuuom  **a|  an 
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Fa  I  A cm*  m»n«  *aaa  *bk  «*ra»  It,  Modi  a*K  o  •  JO  KB "* 
loot  an  tar  ■*■*** Mi  a*u  an  dwwa  a  F|  I 

approximat*  tachtuqu*  with  limited  txpanaon  coaffkamu. 
an  ha*«  ihoaan  ho*  th*  aytmmlwat  and  fimdamantal  mod* 
flaid  fwKtxxa  bahav*  uan|  a  mitt*  component  txpanaon 
and  than  mcraawni  th*  numMi  of  txpanaon  cocffkcKnu  up 
to  four 

Fundamental  mod*  fwM  intcnwry  w*>i  aim  war*  calculated 
for  tiff* ran i  cavity  conflgurttioni  In  particular  aw  haw  fo 
cumd  our  attention  on  th*  lateral  field  ^ot  see*  at  a  function 
of  th*  ectm  rayon  (am  width  Thu  width  can  of  count  br 
aaaocMWd  with  cwnaat  ipreada*  from  th*  nnpa  contact  pro 
tided  th*  cotnpWt*  yometncaJ  p*ct»*»  Fig  1 .  u  known  (1 1 1 . 
1 1 J)  For  *umpM  currant  ytaaduig  from  th*  metallK  ttnp* 
contact  a  *  function  of  th*  ruwtmtm  and  thidumei  of  th* 
*anou>  layran  bale*  th*  conteci  From  Fig  5.  aw  «•  that  th* 
fundamental  mod*  mot  a a*  a  *  function  of  th*  gam  rayon 
width  la,  th*  act m  leyar  thirknaai.  ad  th*  panawtat  ft 
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a* Mar  jam  Ai  i V  (air  m*r«aart  IV  (air  (utdin(  effect  on 
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Ai'-TRAC* 

Injection  li»tr*  having  non)  1n*a-1 1  'es  0* 
'IIMI*  1  r  their  f-I  , POm*'- CuTfn;  Curves  are 
due  primarily  to  lateral  mode  hopping  and  wave¬ 
guide  shifting  with  cu*re«t  This  p«p*>-  H  aimed 
ai  tne  'spates'. to*  of  linea-  devices  wit*  due 
cons ide*at  tons  o*  waveg,ide  aeit?*  fo*  lateral 
mode  stability  **  "CKJr '  a  teterojtrvct or*  laser 
are  calculate  a  lateral  dielectric  ditcontinu' t> 
pairs  o*  tne  gam  loss  profile  o'  tne  setter 
live*  ate  the  overall  geometry.  iMt  ditconti- 
nutty  serves  at  a  'igu*e  of  merit  for  lateral 
aode  stability  we  compe*e  twe  important  double- 
neteroitructure  lasers  (l;  tne  o»'Oe-str'Pe 
geometry.  arc  12  tne  channeled  substrate  plana* 
CS*  .  Calculation*  snow  that  CSf  devices  nave 
strong  'ate*al  inQe*  confinement.  whereas.  tne 
oa'de-stripe  lasers  nave  a  lateral  gat*  confine- 
mtnt  dependent  upon  drive 


;  a ’TOXIC*!  OS 


Tne  ansi.vttt  o'  tne  late*a;  modes  of  stripe 
geometry  laser*  't  complicated  b*  tne  f*ewie*tly 
uncertain  dielectric  0*0*1 'e  in  t*e  Junctto* 
plane  Cu*rent  inte*ett  center*  or  structure! 
tnat  produce  'jndamenta'  lateral  nod*  operation 
ove'  ai  wide  a  pewe*  «m1*si©r  range  at  possible 
we  cwoare  two  type*  of  plana*  Stripe  late'  struc¬ 
ture!  wi  th  the  objective  o*  determining  tne 
structural  requirements  fo*  fundamental  lateral 
node  ooe*ation.  at  analyze  the  oi'de-de'ined 
Stripe  geometry  cpnf 'Juration  (1,2/  and  the 
channeled-sutrstrate  plana*  late*  C 3 ) .  I*  both 
cases.  •*  deal  with  structure*  wne*e  tne  deg-ee 
o*  radiation  tpnead  '*ow  the  active  region  into 
tne  vertical  regions  of  the  device  is  controlled 
Because  there  a*e  regions  o'  varying  loss  tne  P*o- 
pa  gat  to*  constants  for  the  lateral  mdes  dif*e* 

; t  is  possible  to  produce  substantial  dif'e*encet 
betwee*  tne  fundamental  and  tne  nigne-  o*de'  lat¬ 
ere'  nodes  --  a  condition  favorable  fo*  main- 
td'nlng  fundamental  node  operation  unde*  certain 
conditions  discussed  m  this  pape* 

•Supported  1r  part  by  the  U.  S.  »rwy  Research 
Office. 


Tne  theory  is  developed  to  the  eatent  that 
sensitivity  0*  lateral  aodes  to  stripe  width  is 
•  now*  one  a  ar  effective  dielectric  stet  in  tne 
horizontal  direction  is  calculated.  This  di¬ 
electric  stet  H  a  consequence  of  basic  five-lever 
waveguide  considerations  and  is  dependent  up©*  de¬ 
vice  geometry  and  pertinent  materiel  parameters 
Cutoff  conditions  *0*  the  *undanenta1  and  higher 
order  mode*  are  derived  at  a  function  of  the  lat¬ 
eral  dielectric  itet  and  tt'ipe  width.  larller 
wort  related  fundamental  mode  operation  to  wave¬ 
guide  width  only  fo*  tne  strip-loaded  wave¬ 
guide  (4). 

riAVC&UlM  CMRACTtRlSTICS 

*ne  oeta'led  geometry  of  tne  oalde  stripe 
structure  is  Show*  in  *19.  1(a).  regions  denoted 
by  unprimed  n^aoe*s  He  benedV  tne  metallic 
ttnpa  aems  le  p*imed  regions  lie  beneetn  the  oalde. 
*n*  ‘lelds  a*e  confined  vertically  to  region  *  due 
to  tne  'elective  rnae>  barr'e*s  produced  by  the 
two  he te*o.i unctions  at  the  3-4  and  4-S  region 
boundaries  Lateral  confinement  to  the  region  be¬ 
low  the  metallic  contact  is  due  to  the  variations 
in  tne  cample*  d'e'ect'lc  constant,  irfilch  is  re¬ 
lated  tc  tne  gam  and  losses  of  tne  active  liven 
<ae  eit.mw  current  spreading  u-oe*  the  stripe  con¬ 
tact  is  negligible  and  at  a  first  orde*  approai- 
amt  Ion  au/u  the  ga'n  In  region  a  and  tne  absorp¬ 
tion  coe**ldient  In  region  Z  are  constants. 

Since  the  inge*  o'  tne  active  region  is  comt  Heated 
by  anamalous  dispersion,  thermal  and  'nee  Carrie- 
effect* (S;  tne  actua  *e'*act1ve  mae*  ste;  ma»  be 
positive  or  negative  The  geometry  e'  the  IV 
late*  H  tnoav*  m  fig.  1(b)  4nal.nl*  of  this  de¬ 
vice  proceeds  along  lines  similar  to  that  of  the 
oalde  tt'ipe  late*. 

1*  fig  1  tne  dasned  lines  represent  bound¬ 
aries  m  tne  lateral  direction  which  separate  re¬ 
gions  having  different  comrlti  dielectric  con¬ 
stants.  1*  some  waveguidlng  structures  (<  these 
latere'  boundaries  a*e  not  real  but  are  referred 
tc  as  pteude  because  modes  are  confined  when  a 
pseudo- inde>  step  arising  from  different  propaga- 
t lor  constants  on  either  side  of  tne  boundary  be¬ 
comes  positive  The  lateral  boundaries  here  are 
not  pteude  since  the*e  is  a  real  gam/loss  profile 
in  the  active  lave*  that  provides  a  guiding  mech¬ 
anise  even  if  the  lateral  inde«  step  is  negative 
ano  antl-1nde>  guiding  occurs  Asitm'ng  a  mode 
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To  *»C'1lt4t*  ■■Chin*  Calculation!  **  hO»K*11;*  th* 
1*t*r*l  war*  •MK'*-1  and  act  1  »*  -*3'  Or  mitt*. 

"1  *  V»e  *  s  *  V  (•) 

3uto"  *or  •  particular  ’at*ral  km*  occuft  •***- 
4**!(f*4  -  3.  '  *..  th*  KM*  'J  net  gor'ined  tC 

tn*  act«**  -rgion.  uitn  hi  •  n,  .  m*)  tne  setting 
n*  •  C  ■•  »Utn  a  f#t  9'  1n«ol*lnp  S 

ana  in*  ;oaet«i  ai*l*ctric  stec  The  ic  gt'w  of 
dtri  i|«  aquations  ar*  altar  In  FI9.  2. 

>  on) tnal*  14  SA*  1/2  and  the  *64044*  1* 
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act'**  r*c*on  pain,  p,  -row'-rc  'or  tn*  l*t*rt 
at  a  'gnet'or  0'  4  f’p  )  snows  9.»  'or  tn* 

for  C  •  2.C  u»  tn*  9* In  of  tn* 
inc-**t*t  -ap'dly  'o'  1  ’  !  j» 

'or  tn*  treone  nod*  1ncr**t*4  «or 
'■pi  1*4  an  OPt'KK  4  V*lu*  '1*4 
jn  ¥n*r.  o  •  C  5  jt.  gtn  **lu*t 
*-*  n'pn*'  preaut*  H9M  '*a»t  into  tn*  p-6*At  cat 
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for  all  tt-ucturrt  tn*  a*m  cuort  snow  that 
'un )ui»ta i  no at  operation  alwayt  occu-t  oecaut* 
it*  pain  it  tMilttt.  tsourrer.  actual  dnicet 
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*mt  occu-4  oecaut*  pain  pro'M*  dtfonaationt 
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The  Effect  of  Device  Geometry  on  Lateral  Mode 
Content  of  Stripe  Geometry  Lasers 

JOSEPH  B  DLLA\f>  am*  JIROMI  k  BITLER.  siniok  siixihim  it.i 


Tbc  JmIkiik  ptofktr  i.i  ilnpr  puadn  i»|kih«  luer*  a 
(•MxMcd  •ilk  an  ut)Klm  <•!  llIlKlalt  Jnagn  ittiwmifau  lot  lamb 
incaul  Until  muh  ofniliua  Hflnialiaclan  Ua*n  an  aadalnl 
•Ilk  •  MkIik  tup  paortk  HM|  an  tflnlm  Ukvln.  ilaitjalwam 
kairt  oa  Ilka  pia  Uaa  ptufila  o f  Iba  Sell**  la  vat  at  » ail  at  tki  inarall 
piHncinul  ntuciuta  The  laallia  panaadat  a  aaaaiilktnt  naapanNa 
at  iba  patloamaakct  ol  two  aapauai  tkonMr  katanatitklitiA  Uam 
1 1  iba  iiajdn-ttnpa  paokbam  laaat  and  2 1  Iba  cbaiuaUd  anbairau  planar 
l(Vi  laaat  kbodat  iW  uuda  tlnpc  lawn  ba«r  Ulatal  gain  ‘  tmfmamatil. 
whereas  tuklti  ol  CSf  tkevtre*  kiw  kttuakg  lakatal  iwctex  cmIim«<«i 
Ti>  nlatt  Iba  Mkttieebce  ol  gromelrc  t«  Iba  effeclree  Mttint  ptoDIa 
•a  anna  Ibal  Iba  fatal  refract/**  twde*  ol  iba  sen**  la*  at  m  p*aatMie 
tiHlapanttnnl  Kr**iiuik|  tiktUlHaa  u*..»  ibal  i  tlnpa  |tiaaaln  laaat 
tfibataail*  ha»  a  dapnnaad  rffecfi**  inOex  in  INa  act/**  npta  baton 
tba  matallK  contact  Tbit  phenoanekkon  akotaa  ptudwai  inOa t  anti 
pttJinp  la  actual  bncn.  both  paomatr*  aa4  Itta  camat  l*|rcI>on 
iblo  iba  actna  itpoa  pnwlnct  lairtal  lUat  aattfkudiag 
Lalaral  «od«  cat  otl  t of*4t tn no  an  calculated  at  luactuNU  o I  Ik. 
efTactw  complex  dtaWatlK  tup  and  tba  tlnpa  uiallb  Tba  ratull* 
tbou  Ibal  (aloft  a  relate*!  in  a  unique  latbnui  to  Iba  ratio  ol  iba 
taaJ  and  tnaptan  pain  of  iba  complex  Maim  nap  tba  ratio  It 
p«  bailee  foe  in  da  *  pleaded  mode*  and  napaloa  loe  pain  podad  i«t 

Krattn;  ctios 

THf  transverse  mtHies  of  contemporarx  douMe-hetcto- 
function  iatet  diode*  can  be  anal*  fed  rather  accurate)* 
using  a  three  Laver  dielectric  ttructure  In  particular  the  near 
and  far  field*  car  he  calculated  from  the  known  dielectric 
steps  at  the  two  heterotunction  houndartet  and  the  thickness 
of  the  actue  repon  |l)  The  analvxtt  of  the  ferera/ modes  of 
Later  diode*  however  tt  complicated  h*  the  uncertain  die 
eiectnc  profile  in  the  function  plane  o'  ttnpe  (tet'metr*  device* 
Current  interett  center*  on  ttructure*  that  produce  tunda 
menu.  Lateral  mode  operalK«r  »vr-  a*  »ide  a  power  e miction 
range  at  possible 

Previous  analvtet  of  the  latera  mode*  of  tlripe  geometrx 
laaert  ha*e  Been  made  by  assuming  tanoua  dielectric  profile* 
U)-|b|  The  imapnarv  part  of  the  dielectric  constant  was 
tied  to  the  pin  lota  profile  while  the  real  part  wai  related  to 
the  refractive  index  and  it*  associated  perturbations  owing 
to  the  competing  influence  of  ftee  carrier  infection  into  the 
active  laver  and  temperature  profile  |<|  Recent  experi 
mental  measurementt  of  the  lateral  far  field  radiation  pattern* 
of  ttnpe  geometrx  devices  indicate  that  the  local  refractive 
index  of  the  active  repon  ts  depressed  which  impiiei  that  the 
mode  it  strongly  pm  guided  |?| 

Mifiwu  *»>•  rbfbfC  I*er*mbe-  4  jar*  revised  1  ehruar*  2  l®’4 
Tbi*  work  «•<  uipponcd  b»  the  l  lined  Swe*  Arm*  geuir.k  Office 
The  author*  are  ui*h  the  lurifimm-  of  t  metrical  Tnemeerin* 
Southern  Method*' l  m*erun  I bkllu*  T\  *52’* 


In  this  papet  the  lateral  pm  lou  profile  of  the  a* live  laser 
it  modeled  at  a  step  while  lit  tefiactive  index  it  conttam  An 
effective  dielectric  profile  dependent  upon  the  overall  geom¬ 
etry  |S» j  ix  alao  modeled  ai  a  step  Wc  present  a  comparative 
anal* six  of  two  types  ot  planar  stripe  geometry  laser  structures 
with  the  objective  of  determining  the  structure  requirements 
lor  fundaments.  Literal  mode  operation  We  analyse  the 
oxide -defined  stripe  geometry  configuration  and  the  channeled 
lutnirate  planar  tCSf’i  lase*  ft*)  |10]  The  effective  complex  dr 
electru  constant  step  in  the  active  laver  ts  determined  by 
solving  lor  the  propaption  constant*  in  the  transverse  planes 
through  the  contact  and  through  the  oxide  [hie  to  wave 
interactions  in  the  vertical  region*  under  the  contact  and 
oxide  the  effective  dielectric  step  i*  a  strong  function  of  the 
overall  geometr*  at  well  ai  pm  iossei  in  the  active  layer  B) 
assuming  a  uniiorm  index  of  refraction  in  the  active  layet  wt 
isolate  the  influence  of  device  strudure  on  fundamental 
lateral  mode  operation  Calculalrom  show  that  when  the 
active  laver  »  ge omeiricall*  close  to  the  metallic  stupe  con 
tact  the  effective  refractive  index  of  the  active  layer  is  de 
prexted  in  the  region  below  the  contact  thrt  alone  contributes 
lo  index  antiguiding 

The  theorx  is  developed  lo  the  extern  that  dependence  of 
lateral  mode  operation  on  itripe  width  it  known  once  an  eftec- 
uve  dielectric  itep  in  the  active  laver  it  calculated  Thu  ttep 
leadt  to  a  ter  of  cut-off  curve*  for  the  lateral  modes  at  a  func 
iron  of  ttnpe  width  Calculation!  indicate  thar  all  ptn  guided 
lateral  mode*  ate  prevent  under  the  ttnpe  though  then  modal 
gain*  ma\  be  insufficient  fot  oscillation  Lather  work  related 
fundaments  mode  operation  to  waveguide  width  onlx  for  the 
stnp-loaded  waveguide  |ll|  Calculated  far-field  patterns  for 
narrow  ttnpe  iatert  teflect  fundamental  mode  distortions 
desenhed  bv  As  heck  ft  al  |'| 

WAvtc.tinr  CMAHArrrmsTics. 

The  detailed  geometry  of  the  oxide  stripe  structure  ts  shown 
in  Pig  I  in  which  regions  denoted  by  unpnmed  numbers  lie 
beneath  the  metallic  stnpe  while  region*  associated  with  the 
pnmed  numhen  lie  beneath  the  oxide  The  fields  ate  confined 
vertically  lo  repon  4  due  to  the  refractive  index  harriers  pro¬ 
duced  hy  the  two  heterojunctions  at  the  ?-4  and  4-5  repon 
boundaries  Lateral  confinement  lo  the  repon  below  the 
metallic  contact  is  doe  to  the  variation*  tn  the  effective  com¬ 
plex  dielectric  constant  The  mathematical  model  follows 
from  that  first  developed  by  Marcatili  |I2)  and  focuses  alien 
non  on  the  unshaded  repon*  shown  in  Fig  2tst  The  width 
and  thickness  of  repon  4  are  s  and  J.  respectively  The 
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forme:  it  taarr.  nonun*llv  u  ihe  strtpe  width  For  weii-fuided 
mode*  the  field  deefyt  exponential!!  in  reports  3  5  and  * 
The  power  traveling  in  the  i haded  repom  it  negupbK  me  no 
attempt  a  made  lo  match  the  fieidi  aloof  the  edges  of  the 
chided  treat 

Current  flowing  mio  the  acute  laver  frorr  the  contact  spreads 
laieraUs  with  datunce  from  the  contact  ar  that  the  actual 
"fain  width"  may  be  different  from  the  ttnpe  contact  width  t 
The  defree  of  current  ipread  a  determined  by  came:  diffusion 
and  depen dt  on  the  thicknesses  and  rewttvnm  of  the  p-OaAi 
"cap"  layer  and  the  p-AJGaAs  confining  wall  |I3]-|I6|  Fi» 
2(b)  indicatet  the  tier  profile  nf  the«fTecti«e  fam  abaorption 
and  the  effective  in  dee  profile  a  illustrated  m  Fig  2ict  The 
actual  refractive  index  of  the  active  repon  a  complicated  b> 
anomalous  dtsoemon  thermal  aid  free-came:  effects  (51. 
theae  effecti  are  not  included  in  our  inaJssn  since  we  are 
sowing  the  influence  of  overall  device  geometn  on  later 
performance  In  our  model  the  effective  dielectric  profile 
a  modified  rath  hv  wave  interaction  with  the  active  repon 
pm  and  the  overall  ft  onset  rscii  structure  For  example  the 
Aetectnc  profile*  of  the  ttnpe  contact  and  CSP  iaaen  are 
affected  bv  the  metal  oxide  laver  and  the  etched  channels 
respective*! 
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The  optical  field!  are  dinoed  into  two  sets  of  modes  e 
and  E\  ,  |12]  Since  contemporary  heterottructure  lasers 
nave  tlun  active  repons  id)  ~  0  1  to  0.3  iimi  and  electric 
fields  poUrued  predominant!!  aiung  the  lateral  direcuon 
(1*|  l.t  direcuon  in  Fig  2t*i).  we  need  consider  only  the 
E’t  „  modes  Furthermore  we  restrict  discussion  to  even 
mooes  with  res  pec  to  x  ■  0  and  y  ■  0  and  assume  a  mode 
ot  the  form  expo«jr  -  vr  i  where  >  ■  -C  I  •  16  is  the  complex 
propagation  constant  The  £J  s  modes  are  determined  from 
the  field  equations  nhe  superscript  y  u  dropped! 
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wherein  is  the  iree-space  wavenumber,  h,  u  the  free -space 
wave  impedance,  and  a  ts  the  relative  dielectric  constant  £, 
unifies  the  wave  equation 

7}  •  ♦  (*i»  *  7*  id  ■  0  (2) 


where 


Refernng  spin  to  Fig  2(ai,  the  3-d  and  4-5  lines  represent 
teal  vertical  boundaries  separating  two  distinct  materials 
These  boundaries  are  of  course  hetero junctions  with  built-in 
index  steps,  sc  that  light  guiding  results  from  t  positive  index 
increment  The  4-4  dashed  lines  represent  boundaries  ut  the 
lateral  direcuon  which  separate  repons  having  different  com¬ 
plex  dielectnc  constants  In  some  wiveguidmg  structures 
(]R|  these  lateral  boundaries  are  not  real  but  are  referred  to 
as  “pseudo'  because  modes  sre  confined  when  the  effective 
refractive  indei  step  is  derived  from  the  tmapnarv  pans  of 
the  complex  propagation  constants  on  either  tide  of  the 
boundary,  becomes  posltvt  The  lateral  boundaries  of  Fig 
2(a)  separate  gam  and  loas  repon*  m  the  active  liver  so  that  a 
guiding  mechatuun  a  provided  reprdlew  of  a  positive  ot 
negative  effective  index  step  ar  the  boundary 
The  importance  of  the  lateral  4.4  boundary  becomes  *p 
parent  aftet  we  develop  the  etprnequatiofi  pertaining  to  lateral 
confinement  In  repons  4  and  4'  the  solutions  of  the  wave 
equation  become 
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Thr  cottstam  A,  u  determined  from  normalisation  conditions 
lr  the  ah  repon  the  propagation  constants  hJ(  and  h>(  tir 
related  to  the  complex  dielectric  constant  x,  by  (he  equations 

*i.  ♦*},  *  T1  1 4j  i 

hi  -  h\ 5  •  «)•«  ♦  y  i4hi 

Matching  fields  the  4-4  boundary  ate  obtained  bv  placing 
ht'  •  h4>  Thu  pvet  the  lateral  cigcncquauont 


In  1M  the  complex  dielectric  ttep  .1  it  written  at 

A  ■  A’  ♦  iA*  •  »4  •  *«  Ci 

Ihe  imapnarv  part  of  the  dielectri.  constant  of  each  repon 
»  determined  from  elemental  considerations  hoi  example 
with  harmomc  time  variation  expsiwrt  the  simplex  diek.tn. 
conttant  taliet  the  form  »,  >»,  no,  c0wi*x|  n*  where 
n,  a  *J‘  •  u  tt»e  tnJex  of  refraction  o,  u  the  optical  ..undo* 
tivirx  in  the  ah  repon  and  t  tt  the  vacuum  permrltivm  In 
repona  of  optical  absorption  irepor.  lion  positive  while  in 
repont  having  population  tnvetuon  trepor.  4i  o  u  negative 
As  a  remit  the  imapnarv  part  of  ^  »  poutive 

It  u  important  to  note  that  the  dielectric  ttep  in  Ci  it  the 
difference  between  the  dielectric  conitanti  of  repont  4  and  4 
To  include  the  effeett  of  the  light  interaction  of  the  active 
laver  with  all  repont  it  it  ne-risarv  to  define  the  effective 
dielectric  ttep  |10|  |11J  { 1 K J  Thit  lateral  ttep  it  found  bv 
fint  calculating  a  complex  propagation  constant  f  for  a 
transverse  mode  uung  a  five  layer  model  |!)  with  material 
identical  to  that  of  repons  1*5  \nothe:  propagation  con¬ 
ttant  f  b  computed  uung  repont  I  -5  The  effective  complex 
refractive  indexes  of  fhese  modes  are  n  »  i{  and  n  •  -i(’  40 
The  effective  dielectric  step  hecomev  ±  »tf':  f:  i  The 
lateral  mode  character  u  governed  bv  the  effective  ttep  whuh 
includes  the  actual  dielectric  crrutanti  of  the  active  lavet  at 
well  a«  thoae  of  all  vertical  repont 


Cl  TOM  CoNPITIOMS 

The  hauc  condmom  relating  The  waveguide  tmpr  width  to 
the  dielectric  step  at  cutoff  for  the  various  lateral  modes  is 
discussed  in  this  section  For  the  structure  shown  in  Fig  2(a). 
the  eigenvalues  governing  the  field  shapet  of  lateral  modes  are 
determined  from  solutions  of  (5»  for  even  modes  while  the 
odd  modes  satisfy 


h. 


h,  ctn 


(8i 


I  For  simplicity  .  the  v  subscripts  have  been  dropped  » 

In  addition  h,  and  h4  satisfs  (b>  fot  both  even  and  odd 
modes  The  tatio««  *«  a  f  for  contemporary  laser  structures 
To  facilitate  machine  calculations  we  normih/e  the  lateral 
wavenumbers  and  active  repon  width 

(«) 

5  •  k9t  (lOr 


I  i*.  5  <.*ei>e:*iue«l  cutoff  curvei  tor  thr  Utetal  ttrr  mivlcl  ol  an 
imcvtiof  ttte:  Ws>e»u*Sci  with  values  irpreirnicd  pouiti  ir  the 
UiaOrd  irgion  piopofatr  onls  the  lundamenui  lairtaJ  mode 


Cutoff  for  a  particular  lateral  mode  occurs  when  Re  t/f4  I  *•  0. 
i  e  the  mode  is  not  confined  to  the  active  region  Fquationt 
(5  i  and  IM  ate  separated  into  teal  and  imapnarv  parti  uung 
//,  *  H ,  •  t*f  After  minor  manipulations  and  setting  H\  *  0 
we  obtain 
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Fhminating  H\  .  in  the  latt  three  equations  we  have  two 
uniultaneout  equations  with  the  tartiblet  defined  at  H ,S  and 
HlS  where  the  term  £  S’  it  treated  at  a  parameter  In  the 
case  of  lossleu  waveguides  A  A  ■  «  Fot  stripe  geometry 
structures,  the  ratio  £  £  may  of  course  be  either  poutive 
or  negative  Therefore,  for  a  pven  ratio  £  £  the  quantitiei 
H .5  and  H\S  can  be  computed  and  substituted  into  (III  thus 
pvtng  the  value  of  5  at  cutoff  Fig  5  shows  S£‘x  5  as  a  func 
non  of  j’  Waveguides  with  values  represented  by  points 
ir  the  shaded  repon  propagate  only  the  fundamental  lateral 
mode  and  points  above  the  curve  labeled  mode  2  propagate 
the  fundamental  and  second-order  modes 
To  understand  the  meaning  of  these  curves,  it  is  important 
to  understand  the  ortpn  of  the  £  and  A  values  The  real  part 
of  the  dieiectiK  step  is  related  to  the  effective  lateral  index 
step,  it  is  poutive  for  index  guiding  and  negative  fot  gam  ot 
index  antiguiding  The  imaginary  part  it  proportional  to  the 
effective  gam  and  lomes  of  the  active  liver  The  fact  that  the 
ordinate  contains  A  implies  that  fot  o  given  S.  the  gain  of  the 
active  region  can  be  made  large  er.ough  to  make  ant  mode 
propagate  Mode  selection  ts  solely  due  to  the  canty  those 
modes  requiring  the  smallest  gain  from  the  active  region  are 
excited 
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Ir  ’hi*  x.tion  we  diKuts  v>nx  tpecifn  waveguide  geometrxt 
and  calculate  lateral  mode  shapes  and  their  radiation  patternt 
The  ijtffai  m ode  properties  aif  determined  trom  thr  effective 
dielectric  step  at  the  4-4  houndarx  thr  value  of  the  imaginary 
part  of  thr  dxlectru  step  u  dependent  upon  the  gain  lots  of 
regions  4  and  4  at  well  at  thr  overall  tax:  geometry  The 
eftrctive  dielectric  »trp  fot  a  given  lax:  structure  it  uxd  in 
(5l  and  IM  to  obtain  thr  latrra)  eigenvalues  A, ^  and  A(>  for  a 
given  stripe  width  i  and  threshold  gain  ttlh  Thu  active  region 
jam  was  aJtusted  to  given  a  mode  which  travelt  with  G  ■  TO 
cm'1  to  account  for  radiation  can  tv-end  loxxs 
Two  tt ruciuret  examined  here  are  thr  ttnpe  geometry  in 
Fiji  4<ai  and  thr  channel  (CSPi  in  Fig  4ihi  Rrptrxntativr 
matrrul  paramrtrrt  are  included  in  tlx  figures  Of  particuLat 
interrtt  it  thr  ratio  .1  A’  which  it  a  meaturr  ot  thr  strength  of 
index  lot  index  anti-  •  guiding  in  a  lax*  and  u  a  coordinate 
value  for  thr  cut-off  curvet  j  .i’  fot  thr  ttnpe  geometrx  it 
plotted  in  Fig  'tai  versux  d ,  the  p-wall  thickness  Ax  the 
p-wall  shrinks  Oett  than  I  ami  the  defneuting  nature  of  the 
contact  becomes  apparent  in  larger  nrtatnt  valuft  o!  .1  S 
In  thit  model  no  lateral  index  variation  in  the  active  later  it 
attumed  (fain  guiding  only  I  v  that  negative  S  Jt  tatiot  ate 
solely  the  retult  of  geometry  and  proximitx  of  the  active 
region  to  the  contact  The  ratio  £  S'  for  the  channel  ttrux 
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lif  *  I  uim!iuA4  Jcptfttiffivf  oi  thr  dteirctnv  tt(f  tain  Z  Z '  on 
f<  r»r:?>  •  micvttor.  la«m  (it  i  Z"  %rtwi«  the  confining 

•hhinrtt  It*:  *  fcenwtn  la»cr  «n0  (hi  ^  z  %crui»  thr  tor 

■  inm*  r  mai  thKlnen  lot  »  (  SI*  i iw? 

turc  of  I  ig  4<pi  ii  piv-tted  ir,  Fig  Shi  .ertut  the  confining 
n  wall  thickness  T tie  -hanne'  depth  it  1  T  um  til  patametett 
ate  ..hoxr  to  cloxlx  rrxmhie  thox  of  Aiki  e;  el  |10|  The 
tatii  j  ^  ii  potitivr  when  thr  guiding  nature  of  ihr  channel 
it  enhanced  for  f  <0  .5  um  The  channel  mode!  of  Fig  4<bi 
d.xt  not  include  the  influence  of  the  contact 
Active  regior  gain  flf  it  plotted  at  a  function  of  ttnpe  width 
j  for  the  ttnpe  geometry  and  C SP  model)  in  Figs  Na i  and 
•  hi  respectively  In  Fig  oral  loi  d,  *  I  0  um  the  gain  of 
the  tundamema  mode  meteaxt  rapidlx  for  i  <  *  i/rr.  whereat 
the  gait  tot  the  xennd  mode  metraxt  for  t  <  10  um  Thit 
implies  an  optimum  i  valur  Ixt  between  5  and  10  *im  When 
cf,  •  0  5  arm  f„  values  are  higher  hecaux  light  leaks  into  the 
p-OaAt  cap  In  Fig  whi.  gltl  fot  the  fundamemal  mode  of 
the  CSP  laxr  hat  a  tofi  dependence  on  /  between  0.T  and 
OF  urn  However,  fot  the  xcond  mode  glh  nxt  rapidly  fot 
»  <  10  am i  and  r  •  0  F  am  For  favorable  tingle  mode  opera 
tion  we  conclude  that  fot  r  *  0  J  aim.  t  <  “  am  while  fot 
I  •  O  F  jam.  »  <  F  um 

Fot  all  structures  the  gain  curvet  show  that  fundamental 
mode  operation  always  occurt  hecaux  itt  threshold  gain  it 
smallest  However,  actual  devtcet  with  large  i  salues  opetate  in 
high -order  modes  Thu  occurs  hecaux  gain  ptofile  deformi 
tiont  caux  switchings  of  mode  preference  Thex  instabilities 
are  frrquentlv  obxrved  in  ttnpe  laxrt  and  to  a  letxr  degree  in 
CSP  devicet  The  ptomnxnt  reason  fot  trrstabilitxs  it  the  fact 
that  ttnpe  geometn  laxtt  have  gain  guided  modet  influenced 
hx  gam  profile  delormationt  Thex  index  antiguidrd  ms»det 
haxe  negative  S  S'  valuex  Fven  though  the  actual  rcltactixr 
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index  of  the  active  layer  »  uiumtc  constant  the  rfftcmt  8<bi  where  d,  •  0.?  um  the  twm  peak*  are  more  width 

index  uep  u  negative  becau*  of  mode  interaction*  with  the  ‘paced  The  iarger  negative  £i£"  value  for  d,  •  0-5  urr.  ha* 

metallic  and  oxide  conucta  a  defocusng  effect  on  the  fa/ -field  pattern  The  twin  peaki 

For  a  giver,  itnpe  width  a.  the  values  of  5A*,/:  and  A  £"  are  in  these  pattern  are  urruiar  to  partemi  calculated  by  Aabeck 

calculated  Contoun  repeesenung  the  different  atructuret  are  era/  (7)  for  narrow  nnpet  uatng  an  Eckan  potential  model 

computed  lot  the  fundamental  mode  gain  In  Ftg  7  w»  plot  Far  field  pattema  lot  a  CSP  laaer  are  ahown  in  Fig  8tc)  and 
the  locua  of  ootnti  on  the  cutoff  curvet  High- order  mode  <di  In  the  limit  ei  the  n-confinmg  wall  ttuckneu  t  get t  large, 

contour*  would  follow  those  of  the  fundamental  ai  large  £  A*-*0  Hence,  the  pattern*  of  Fig  8(c)  ere  armlet  to  those 

or  dm  ant  value*  but  a*  r  *•  0  the  contoun  approach  different  of  Fig  Mai  where  A  /A*  u  small  In  Fig  8(d).  /  ■  0J  um  and 

limiting  point*,  foe  example  the  second  mode  terminate*  on  the  twin  peak*  are  absent,  the  light  ta  closely  confined  to  the 

S£*'  ’  •  3  ’5  and  £  A*  ■  0  Note  that  for  the  nnpe  geom  center  of  the  nnpe  where  the  defocutmg  effect*  of  the  light 

etrv  device  with  d,  •  0  5  um  t curve  hi.  the  fundamental  mode  propagating  in  the  pruned  region*  are  absent  The  focusing 

contour  crosMt  the  *econd -order  mode  cutoff  tt  «  =»  S  ion  nature  of  the  channel  ta  reflected  ut  a  large  positive  A,  A” 
below  thu  r  value  the  second  mode  does  not  exist  at  funda-  ratio 
mental  mode  threshold  gam*  On  the  other  hand,  the  CSP 

with  r»0J  arm  (curve  d>  crosses  the  second  mode  cutoff  at  Discussion 

i  a*  2  J  um  The  theoretical  analytu  presented  in  this  paper  shows  the 

The  £,£"  figure  of  merit  for  vinous  structures  u  demon  powerful  effect  on  the  lateral  mode  content  of  lasers  due  to 

strable  in  farfleld  calculations  in  the  lateral  plane  where  the  the  radiation  extending  substantially  beyond  the  re  com  bin* 
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uon  region  either  towarc  the  device  surface  01  substrate 
Earlier  work  [4|  snowed  the  effect  on  the  threshold  current 
density  ttorr  absorption  ir.  the  external  region i  ol  double 
rteterojuncuon  lasers  The  preient  analysts  indicate!  the  sttong 
effect  on  laterai  mode  propagation  b\  tunable  control  of  ttie 
perpendicular  device  profile  o:  neteroiunction  structures 

The  oxide  -defined  ttnpe  |eometrv  laser  u  theoretical!* 
shown  to  be  sensitive  to  tne  traction  of  the  radiation  rea.lung 
the  tur'.ace  region*  of  the  device  with  regard  to  the  prooaga 
uor  constant*  ot  the  vanout  lateral  modes  The  effect  on  the 
dielectric  constant  u  such  that  anttguiding  of  the  .fundamental 
mode  occurs  Because  of  tne  sensitivity  of  the  propagation 
constants  to  the  extent  of  radiation  spread  toward  the  surface 
of  the  device  land  lack  of  stabilization  of  the  fundamental 
mode  i  changes  in  laser  current  will  produce  perturbations  in 
the  dielectnc  profile  Owing  to  the  tree  carrier  density  dis¬ 
tribution  associated  gam  profile  and  local  heating,  the  cutoff 
condition!  for  various  orde-  lateral  modes  will  change  with 
drive  As  described  m  kirxbs  ei  j i  (4  J  lateral  mode  shifts 
with  current  are  indeed  commonly  icen  in  planar  itnpe  lasers 
Variations  among  lasers  mas  be  partialis  explained  by  dif¬ 
ferences  in  the  degree  of  radiation  confinement  m  the  direc¬ 
tion  perpendicular  to  the  junction  plane 

The  radiation  fteldi  calculated  show  the  influence  of  geom¬ 
etry  on  the  pattern  shapes  In  the  stnpe  geometry  structure 
pattern  distortion  occurs  for  narrow  stripes,  this  distortion 
reflects  the  amount  of  wave  energy  outside  the  active  region 
In  tne  strong  index  guided  CSf  4scrs.  pattern  diatorlion  is 
absrnt  It  should  be  noted  that  these  calcinations  were  madr 
vnlh  a  constant  real  refractive  index  us  the  active  layer  and 
that  index  anuguidmg  results  trom  the  effective  dielectnc 
profile  Inclusion  ol  s  real  negative  index  step  accentuate! 
the  distortion  and  spreads  the  pattern 

The  dielectric  ratio  u  u  positive  for  CSf*  lasers  while  it 
is  negative  for  itnpe  geometrx  devices  CSf  waveguidei  pro 
duce  lateral  confinement  in  the  absence  of  gain  losses  For  a 
reactive  CSP  waveguide  the  complex  propapnon  constant 
is  imaginary  In  our  model  this  results  tn  Imaginarvi  *  I  > 
Imaginary  f '  >  with  the  difference  being  the  product  of  the 
tree  ipace  wavenumber  and  the  effect.'**  index  step  If  the 
.hannei  were  inverted  Imaginary  f  i  <  Imagsnarvif  and  the 
effective  index  step  would  be  negative  there  would  be  no 
lateral  mode  confinement  X  stnpe  geometrx  waveguide  with 
<  perfect  conducting  stnpe  would  behave  similarly  to  an  in  - 
•tried  channel  waveguide  ut  the  sense  that  lateral  mode  con¬ 
finement  would  not  exist 

F  mails  the  step  profile  approximation  loses  accuracs  as 


the  stnpe  width  narrows  Consequently  Figs  b  ~  and  t 
should  be  interpreted  on  a  broad  qualitative  basis  lor  small 
i  values 
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ABSTRACT 

Dispersion  characteristics  and  radiation  fields  of  an  optical  strip- 
line  waveguide  radiating  into  free  space  are  calculated.  The  waveguides 
are  fabricated  as  multiple  layers  of  differing  dielectric  materials.  A 
top  layer  is  etched  to  form  a  "cap”  with  an  effective  waveguide  in  a  lay¬ 
er  below  the  cap.  Confinement  of  the  fields  to  the  waveguide  is  obtained 
in  the  vertical  direction  by  dielectric  discontinuties ,  while  lateral  con¬ 
finement  occurs  because  of  spatial  interference  of  a  continuum  of  plane 
waves.  The  radiation  field  of  the  fundamental  mode  in  a  plane  perpendicu¬ 
lar  to  the  waveguide  levers  is  characteriied  by  the  laver  widths  and  index 
discontinuities.  Beaawidths  of  the  fundamental  node  in  the  plane  parallel 
to  the  dielectric  lavers  are  developed  in  terms  of  the  waveguide  parameters 
Values  of  these  parameters  which  yield  the  best  optical  confinement  under 
the  stripe  can  be  obtained. 


*  Supported  in  part  by  the  U.S.  Annv  Research  Office 
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I.  Introduction 


The  optical  stripline  waveguide  has  potential  applications  as  a  low  loss 
channel  waveguide  in  integrated  optical  circuits  [1].  The  interest  in  this 
waveguide  is  due  to  the  fact  that  it  is  relatively  staple  to  fabricate  cooper¬ 
ed  to  burled  waveguide  structures  having  built-in  dielectric  steps  in  both 
transverse  directions.  Tne  basic  structure  of  the  optical  stripline  is  shown 
in  Fig.  1.  A  similar  structure,  with  region  4  replaced  by  a  conducting  plane 
and  d-“C,  has  beer,  proposed  for  use  at  oilllaeter  and  aubcilllaeter  wavelengths 
[2].  Current  research  is  aimed  at  producing  circuit  elenents  and  system  appli¬ 
cations  using  these  waveguides. 

An  overview  of  dielectric  waveguides  for  oicrowave  Integrated  circuits  has 
beer,  given  by  Knox  [3).  Among  the  aost  recent  applications  are  scannable  anten¬ 
nas  and  tunable  filters.  Structures  have  been  fabricated  for  use  as  electronic 
phase  shifters  at  oilllaeter  [4]  and  subaillioeter  waveienghts  [5].  Itoh  and 
Hebert  .6]  have  sloulated  an  elect ronlcal! v  scannable  antenna  structure  with  a 
mechanical  scan.  A  review  of  work  in  integrated  optics  has  been  given  by  Kogel- 
nik  oanv  of  the  devices  demonstrated  use  the  optical  stripline  waveguide 

as  a  transmission  element.  An  example  is  the  directional  coupler  18),  which 
consists  of  several  striplines  spaced  closely  together  to  allow  coupling  between 
guides  which  has  beer,  analvzed  using  coupled  mode  theories  [9). 

The  effective  index  method  )lj  is  a  mathematically  simple  way  of  predicting 
some  of  the  properties  of  optical  strlpllnea.  This  method  explains  light  confine¬ 
ment  under  the  stripe  by  assuming  an  effective  index  of  refraction  step  in  the  lat¬ 
eral  direction  in  region  3  of  Fig.  1.  This  effective  index  step  is  found  by  first 
obtaining  the  propagation  constants  for  the  structure  shown  in  Fig.  1  with  w  •  * 
and  then  with  w  •  0.  The  difference  between  the  two  propagation  constants  divided 
bv  the  free  space  wavenumber  gives  the  effective  index  step. 


The  far  field  radiation  pattern  of  a  cleaved  striplin*  radiating  into  free 
apace  la  a  useful  wav  to  characterize  the  atructure.  For  example,  the  pattern  can 
he  measured  exper isentally  and  cotapared  with  these  computed  reaults  to  obtain 
values  of  the  various  waveguide  parameters  such  as  the  refractive  indices  n, .nj. 
n,  and  n  ,  or  the  dimensions  d,,  d -  and  w.  Another  potential  application  cf  the 

J  H  *  J 

atrlpline  structure  la  In  optical  arravs  where  the  waveguide  la  used  as  a  member 
of  the  array.  The  pattern  of  the  arrav  of  strlpllnes  la  the  product  of  the  pat¬ 
tern  of  a  single  element  and  the  arrav  pattern. 

In  this  paper  we  develop  a  series  of  rigorous  calculations  showing  how  the 
far-fieid  radiation  pattern  behaves  as  a  function  of  various  waveguide  parameters. 
In  particular,  we  calculate  the  halfpower  beamvidth  In  the  lateral  direction  of 
a  waveguide  mode  radiating  Into  free  space.  The  value  of  the  beaawidth  is  close- 
lv  related  to  the  effective  lateral  index  step  and  conseouent ly ,  determines  the 
degree  with  which  the  waveguide  mode  is  confined  to  the  region  below  the  stripe. 

For  example,  broad  patterns  are  due  to  stronglv  confined  modes.  (These  modes 
should  be  less  vulnerable  to  losses  introduced  by  lateral  waveguide  bends). 

A  rigorous  mathematical  model  of  the  optical  stripline  [10]  does  not  require 
use  of  an  effective  lateral  index  step.  An  extentlon  of  this  model  has  been  de¬ 
veloped  for  the  structure  in  Fig.  1.  The  height  is  assumed  to  be  Infinite.  Num¬ 
erical  results  are  given  to  show  how  the  pattern  depends  on  the  parameters  d^.d^.w, 
and  the  index  discontinuity . 

Since  the  effective  index  approximation  Involves  considerably  less  computa¬ 
tion.  a  method  for  obtaining  information  on  the  radiation  pattern  In  the  lateral 
direction  using  this  approximation  would  be  useful.  One  possible  approach  la  to 
compute  the  effective  Index  step  in  the  lateral  direction  and  then  compute  the 
radiation  pattern  of  a  symmetric  three  laver  waveguide  srtth  this  Index  step  and 
a  thickness  equal  to  the  stripe  width.  We  have  used  this  approximate  technique 


II.  Theory 

The  analyst*  of  the  waveguide  structure  will  parallel  that  given  In 
Ref.  [1C].  Ue  will  assuae,  following  Marcatilli  Ill},  that  there  are  two 
set*  of  nodes,  one  polarized  along  the  x  direction  and  one  along  the  y 
direction.  For  sinpllclty,  we  will  deal  with  the  nodes  polarized  along  the 
x  direction.  Ue  will  further  restrict  this  discussion  to  include  only  the 
ever,  nodes.  The  general  solutions  tc  Che  wav*  equation  in  the  various  re¬ 
gions  can  he  written: 


»'.(«,?)•  _  A  co*  p'  x  ex;  (-p  "(y  -  d,)> 

1  n  n  n  2 


•  •.(x.y)  -  cos  qx  !E(q>  exp(s"(v-d,))  +  F(q)  exp(-s”(y-d,)) )  dq  (lb) 

*  4  4 


v3(x.v)  • 


»4(x,y)  - 


cos  qx  ,'B(q)  sin  q''v  ♦  C(q)  co*  q"r)  dq 

J 

0. 

cos  qx  lD(q)  ex;  r"  ♦  >•))}  dq 


where 


i,  ■  solution  of  the  wave  equation  in  the  1  region 

k^  •  free  space  wavenunber 

n,  •  index  of  refraction  in  the  ilh  region 

k  •  k  n 
1  o  i 


pn  "  2v 


£  ■  propagation  constant 
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Matching  the  field*  and  their  derlvat 
equation  of  the  fora: 


(4a) 

(4b) 

(4c) 

(4d) 

•  at  the  boundaries  yields  a  secular 


Det  (E  -  D  -  0  (5) 

where  E  is  the  unit  matrix  and  T  Is  a  matrix  whose  elements  depend  on  the  un¬ 
known  propagation  constant  8.  Since  the  matrices  In  (5)  are  of  Infinite  order 
the  eigenvalue  solutions  are  found  by  taking  an  increasing  matrix  order  be¬ 
ginning  with  a  first  order  matrix.  This  process  forms  a  convergent  sequence, 
converging  to  the  true  value  of  i .  For  the  numerical  results  given  In  this 
paper,  a  fourth  order  matrix  was  used.  This  order  Is  sufficient  since  the  se¬ 
quence  converges  rapldlv. 

The  constants  B.C.D.E,  and  F  of  Eq .  (1)  are  found  by  matching  the  fields 
or  derivatives  at  the  boundaries.  The  field  description  can  then  be  found  by 
substituting  these  constants  Into  Lq.  (1),  giving: 

»,<x.y)  •  '  A  cos  p’x  exp(-p"  (v-d,))  (6) 

a  n  n  n  i 

n»l 

».(x,y)  ■  2  A  I  i  -  2.3,4  (7) 

1  nil  °  1 

where  represents  an  Integral  expression. 

The  radiation  pattern  Is  related  to  the  Fourier  transform  of  the  field  at 
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the  aperture  112).  Denote  this  transfonr.  by: 

•(Wx,ky)  -  E,(x,y)  e:(kxx  *  kyy)  dxdy 


(8) 


where  r,  represent*  the  transverse  facet  electric  field,  and 


k  •  k  cos  d  sin  6 
x  o 


k  •  k  sin  is  in  6 

v  o 


(9a) 

<9b) 


where  t  denote*  the  angle  f  roe  the  z-axls  and  C  i*  the  angle  froc  the  x~axis 
in  the  x-v  plane.  The  far  field  pattern  for  the  case  of  an  aperture  field 
polarised  in  the  x  direction  can  be  written: 


E(6,«)  ■  C[6*x  cosC  «  de^  sin$  cost] 


(10) 


where  C  is  a  constant  and  6  and  d  are  unit  vectors.  It  it  therefore  necessary 
to  obtain  e  (k  , k  )  which  is  the  Fourier  transform  of  the  field  distribution 

XX? 

given  previously.  The  result  can  be  written: 


,x(kx*kv)  •  I  r(*i} 


(11) 


where  F(«i>  denotes  the  Fourier  transform  of 
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1.  Numerical  Result* 


7h«  propagation  constant  for  the  structure  In  Fig.  (1)  with  d.  •  0  has 
beer,  previously  calculated  [10].  The  effect  of  region  2  on  the  propagation 
constant  is  shown  in  Fig.  2.  We  have  used  normalized  value*  of  the  paraaeters 
so  that  the  plot  will  be  applicable  to  a  wide  class  of  structures.  The 
normalized  paraaeters  arc: 


2 

,  1/2 

"  d:ko(n3 

-  > 

1  /? 

* 

•>  1  '  * 

D3- 

d3ko(n3i  - 

2 

2  1/2 

w  - 

W  fc0(n3*  - 

v> 

2 

.2 

n 

B* 

H 

2 

*0  (n3*  * 

\  )  n3*  -  nA 

(12a) 

(12b) 

(12c) 

(12d) 


(12*) 


(12f ) 


Fig.  3  show*  the  normalized  propagation  constant  as  a  function  of  2V, 
the  normalized  stripe  width,  for  various  values  of  D,.  The  propagation  constant 

4 

approaches  that  of  a  three  layer  waveguide  for  sufficiently  large  values  of 
the  strip*  width.  For  smaller  values  of  the  normalized  stripe  width,  which 
are  more  typical  of  the  values  used  in  integrated  optics,  the  propagation 
constant  is  a  strong  function  of  Dj*  For  large  enough  values  of  Dj,  we  would 
simply  obtain  a  three  layer  waveguide,  and  the  effect  of  the  stripe  would  be 
inconsequential. 
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Ue  now  turn  our  attention  to  the  radiation  fields.  The  pattern  of  a  wave¬ 
guide  radiating  into  free  space  depends  on  the  actual  values  of  the  various 
waveguide  parameters  and  thus  cannot  be  expressed  in  terms  of  normalised  para¬ 
meters.  Consequently,  we  chose  values  of  the  various  parameters  typical  of 

those  used  in  integrated  optics.  Consider  a  GaAs  -  AlGaAs  waveguide  with 

nl  "  n2  *  n4  *  3,55 
n3  -  3.6 

d^«  1.0  um 

Ue  have  chosen  a  wavelength  of  1.15  bm  corresponding  to  one  of  the  transitions 
lr.  a  HeNe  laser.  At  X  •  1.15  us  absorption  in  GaAs  is  low. 

The  halfpower  beaawidth  of  the  fundamental  mode  in  the  lateral  direction 
is  plotted  as  a  function  of  the  stripe  width  in  Tig.  3.  A  scan  in  the  lateral 
direction  was  chosen  because  the  far  field  beamwidth  will  indicate  how  well 
the  fields  are  confined  laterally.  This  is  significant  because  confinement 
of  the  fields  lr.  the  lateral  direction  does  not  occur  because  of  an  index 
discontinuity,  and  sc  it  is  important  to  understand  the  waveguide  parameters 
which  yield  the  greatest  lateral  optical  confinement.  The  degree  of  optical 
confinement  is  Important  in  considerations  of  scattering  at  waveguide  bends. 

From  Fig.  3,  it  is  apparent  that  the  beamwidth  decreases  as  dj  Increases. 
For  small  values  of  the  stripe  width,  the  beamwidth  is  a  strong  function  of 
d^.  Each  plot  exhibits  a  maximum  in  the  beamwidth,  which  would  correspond 
to  an  optimum  value  of  w.  Input  coupling  to  the  waveguide,  which  depends  on 
the  numerical  aperture,  would  be  maximised  for  this  value  of  w.  For  dj  ■  0, 
for  example,  this  optimum  value  would  be  a  stripe  width  of  approximately  3b. 


9 


In  obtaining  Fig.  4,  a  structure  was  used  with  n,  •  a,  •  n^,  tv.  •  3.6, 
d  ■  lb,  X  •  1.15b,  v  •  3b,  and  the  difference  n.  -  n,  was  denoted  &n.  This 
plot  shows  that  the  dependence  of  the  lateral  beaawldth  on  the  transverse 
In  is  enhanced  by  the  presence  of  region  4.  Since  a  large  &n  is  desired  for 
confinement  in  the  vertical  direction,  this  is  the  region  of  the  plot  which  is 
of  Interest.  Note  that  when  in  •  0.1,  the  effect  of  the  thickness  of  region  4 
increasing  from  0  to  0. 2„  is  to  decrease  the  halfpo,*»r  beaawldth  by  about  35*. 

Figure  5  shows  the  dependence  of  the  lateral  beaawldth  on  the  thickness 
of  the  waveguide  region  2.  The  structure  has  n.  ■  a,  •  n .  •  3.55,  n.  •  3.6, 

X  •  1.15,  and  w  •  3b. 

It  is  of  interest  to  note  that  the  effective  index  step  approximation  dis¬ 
cussed  earlier  yields  good  results  for  the  lateral  half  power  beaawldth.  Fig¬ 
ure  6  shows  plots  of  half  power  beaawldth  versus  d.wlth  d„  »  0  computed  using 
the  theory  above  (curve  1)  and  using  the  effective  index  approximation  (curve  2). 
It  is  apparent  that  the  approximation  predicts  the  correct  shape  of  the  curve. 

The  error  which  can  be  made  using  the  approximation  is  less  than  26*.  If  the 
pattern  is  needed  only  in  the  lateral  direction,  the  use  of  the  approximation 
mav  be  justified  by  the  savings  in  computation  time. 
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IV.  CONCLUSION 


The  propagation  constant  t  of  tn«  fundamental  node  in  the  optical  strip- 
line  waveguide  has  been  calculated.  For  snail  stripe  widths  the  value  of 
s  is  strongly  dependent  on  the  presence  of  a  region  between  the  stripe  and  the 
wavegulding  region.  As  the  thickness  of  the  isolating  layer  increases,  the 
effect  of  the  strip  width  on  i  decreases. 

The  halfpower  beasvidrh  of  the  fundamental  node  of  the  optical  stripline 
waveguide  in  a  plane  parallel  to  the  waveguide  layers  has  been  obtained;  this 
is  a  useful  indicator  of  the  degree  of  optical  confinement  in  the  lateral  direc¬ 
tion.  As  the  stripe  widens  the  beanwidth  decreases;  Indicating  that  the  fields 
are  spreading  in  the  waveguide.  As  the  stripe  width  is  decreased  beyond  a  cri¬ 
tical  point,  the  beanwidth  begins  to  decrease;  this  indicates  that  the  optical 
flleds  are  no  longer  confined  uder  the  stripe.  The  effect  of  increasing  the 
thickness  of  region  2  is  al wav*  to  reduce  the  halfpower  beanwidth. 

The  lateral  beanwidth  also  depends  on  the  thickness  of  the  wavegulding  region, 
iy  As  d.  is  decreased  the  beanwidth  Increases,  up  to  the  point  where  d^  becomes 
so  snail  that  the  fields  begin  to  spread  outside  the  wavegulding  region.  Then 
the  beanwidth  drops  sharply.  This  behavior  is  also  predicted  by  using  the  effec¬ 
tive  index  approximation  to  determine  an  index  discontinuity,  and  then  calculating 
the  far  field  beanwidth  of  a  three  layer  waveguide  with  this  index  step  and  a 
thickness  of  2w. 
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FIGURES 


Fig.  1 
Fig.  2 

Fig.  3 

Fig.  4 

Fig.  5 

Fig.  6 


Croat  section  of  the  optics!  stripline  waveguide 
Normalized  propagation  constant  of  the  fundanental  node  as  a 
function  of  normalized  stripe  width 

lateral  halfpower  beanwidth  of  the  fundanental  mode  as  a  function 
of  stripe  width 

lateral  halfpower  beanwidth  of  the  fundamental  mode  as  a  function 
of  the  index  discontinuity 

lateral  halfpower  beanwidth  of  the  fundamental  mode  versus  the 
thickness  of  region  3,  the  region  where  light  is  confined. 
Conparison  of  the  method  given  in  this  paper  for  calculating  the 
lateral  halfpower  beanwidth  (curve  1)  with  the  effective  index 
method  (curve  2) 
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Evaluation  of  Dielectric  Optical 
Waveguides  From  Their  Far  Field  Radiation  Patterns* 
by  M.  V.  Scott  and  J.  K.  Butler 
Southern  Methodist  University 
Dallas,  Texas  73205 

abstract 

At.  experimental  technique  for  determining  the  characteristics  of 
dielectric  optical  waveguides  is  presented.  The  far  field  pattern  of  a 
passive  waveguide  excited  by  an  external  source  is  aeasured  and  compared 
with  a  theoretical  pattern.  The  presence  of  spurious  light  in  the  far 
field  complicates  the  measurement ,  but  this  difficulty  car.  be  minimized 
for  aanv  structures.  Application  to  semiconductor  laser  structures  is 
discussed.  Determining  the  wavegulding  properties  of  Isser  structures 
prior  to  processing  car  result  in  the  early  detection  of  faults  for  in¬ 
creased  yield. 
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Evaluation  of  Dielectric  Optical 
Waveguide*  Froe  Their  Far  Field  Radiation  Pattern* 

Dielectric  optical  waveguide*  are  the  basic  building  block*  used  to  construct 
integrated  optical  circuits  [1,2].  Semiconductor  laser*  consist  essentially  of 
dielectric  optical  waveguides  with  optical  ’’gain"  ir.  the  active  region  to  furnish 
light  amplification  [  >-  7 ) ,  The  properties  of  dielectric  waveguides  depend  on  the 
geoaetry  of  the  structure  and  the  index  of  refraction  discontinuities  between  the 
region  where  light  is  confined  and  the  surrounding  regions.  Device  geometrical 
dimensions  car.  be  studied  using  the  scanning  electron  microscope  (SEX)  and  the 
index  of  refraction  profile  can  be  estimated  from  published  data  [8],  [9]. 

In  this  paper  we  present  an  experimental  technique  for  determining  the  pro¬ 
perties  of  certain  optical  waveguide  structures.  This  technique  can  be  applied 
to  semiconductor  laser  wafers  before  thev  are  cleaved  tc  form  individual  delvces. 
The  wevegulding  properties  of  the  structure  can  be  determined  prior  to  actually 
uklng  the  device.  The  detection  of  any  faults  in  the  wafer  at  this  point  would 
result  in  a.:  increased  yield. 

The  experimental  apparatus  used  to  excite  a  propagating  mode  in  the  wave¬ 
guide  has  been  used  previously  [10-1“',.  Light  from  a  1.13  vn  HeSe  laser  is 
focused  into  the  waveguide  by  a  microscope  objective.  A  second  microscope  ob¬ 
jective  Is  focused  on  the  output  side  of  the  waveguide.  A  camera  with  an  infra¬ 
red  vldlcon  ie  then  used  to  observe  the  near  field  of  the  waveguide. 

After  a  propagating  mode  has  beer,  excited,  a  germanium  detector  with  a  small 
active  area  is  scanned  across  the  far  field.  An  example  of  e  far  field  pattern 
measured  experimentally  is  shown  vith  solid  lines  in  Fig.  1.  This  pattern  was 
obtained  from  a  multi-laver  structure.  An  n-type  lever  of  AlGaAs  doped  with 
tellurium  is  grown  on  a  GaAs  substrate.  This  layer  has  about  11  aluminum;  the 
next  layer  is  GaAs  lightly  doped  with  silicon.  The  final  layer  ie  p-type  AlGaAs 


2 


doped  with  germanium  with  about  201  aluminum.  Tne  thicknesses  of  the  various 
regions  are  known  quite  accurately  from  measurements  with  the  SD:. 

The  middle  layer,  which  it  light ly  doped  GaAs,  will  confine  most  of  the 
light  since  it  is  the  region  with  the  highest  index  of  refraction.  The  refrac¬ 
tive  index  of  this  region  can  be  extrapolated  from  data  published  previously 
[6].  If  this  value  is  not  accurate,  it  makes  very  little  difference  since  the 
waveguldlng  properties  of  the  structure  are  determined  by  the  index  steps  be¬ 
tween  this  region  and  the  other  lavers . 

One  other  observation  can  further  simplify  the  problem  of  obtaining  a 
theoretical  pattern  tc  match  the  one  recorded  experimentally.  Since  there  is  a 
large  index  step  between  the  active  region  and  the  layer  with  201  aluminum, 
there  will  be  very  little  light  confined  to  this  201  aluminum  layer.  We  can 
therefore  estimate  the  step  between  the  active  region  and  the  20 1  aluminum  layer 
to  be  0.1  using  the  references  cited  earlier  (6,9]. 

The  index  of  refraction  profile  resulting  from  the  estimates  given  above 
is  shown  in  Fig.  2.  The  thicknesses  of  the  various  regions  are  also  shown.  The 
index  step  between  regions  2  and  3  and  regions  1  and  3  will  have  a  significant 
effect  or.  the  far  field  pattern.  These  index  steps  are  labeled  An,  and  An^ ,  re¬ 
spectively,  in  Fig.  2. 

For  giver,  values  of  An,  and  Ar...  the  far  field  pattern  of  e  multi-layer 
wevegulde  with  the  index  profile  shown  in  Fig.  2  can  be  calculated  [15], 

Fig.  3  shows  how  the  helfpower  beaavldth  of  the  fundamental  mode  of  the  structure 
in  Fig.  2  varies  es  s  function  of  An.,  for  certain  values  of  An,. 

The  experimental  pattern  discussed  earlier  is  obviously  the  pat  tern  of  a  fun¬ 
damental  mode,  since  it  hae  only  one  peak.  The  half power  beasnfldth  of  the  exper¬ 
imental  pattern  is  1?*.  From  Fig.  3  we  note  that  there  are  an  infinite  number 
of  combinations  of  An.  and  An,  which  would  produce  the  same  halfpower  beamwidth. 


Fortunately,  all  of  these  parameter  values  do  not  produce  the  sane  radiation  pat¬ 
tern.  For  example,  two  radiation  patterns  with  the  sane  halfpower  beamvidth  but 
different  shapes  are  shown  superimposed  on  the  experimental  pattern  In  Fig.  1. 

One  of  these  patterns  agrees  with  the  experimental  result  only  at  the  halfpower 
points,  with  no  agreement  at  other  points.  The  second  pattern  more  closely  matches 

the  experimental  pattern.  The  values  of  An,  and  An.  resulting  in  this  pattern 

1  i 

were  found  by  requiring  the  theory  to  match  experiment  at  the  half power  and  70S 
power  points.  Since  small  perturbations  In  the  parameters  An,  and  An,  result  In 

significant  changes  lr.  the  pattern,  the  parameters  An,  and  An.,  have  been  determined 

1  «* 

quite  accurately  by  this  procedure. 

The  procedure  outlined  above  does  have  certain  limitations.  It  works  best  on 
structures  designed  for  single  mode  operation,  l.e.,  small  active  regions  and  small 
index  steps,  An^.  Although  the  method  car.  be  applied  to  multimoded  structures,  ob¬ 
taining  a  theoretical  natch  to  the  experimental  pattern  is  often  complicated  by  the 
presence  of  several  modes  adding  together  to  form  the  aperture  field.  The  fields 
of  each  of  these  modes  must  be  multiplied  by  complex  weighting  coefficients  to  ac¬ 
count  for  phase  differences  between  the  modes.  This  introduces  several  more  un¬ 
knowns  into  the  problem,  and  makes  the  "trial  and  error"  approach  to  pattern  match¬ 
ing  verv  tedious. 

One  other  potential  problem  Is  the  presence  of  scattered  light  in  the  far 
field.  This  problem  is  minimised  by  using  an  input  microscope  objective  with  a 
spot  size  smaller  than  the  thickness  of  the  active  region  of  the  waveguide.  This 
ensures  that  most  of  the  light  will  be  focused  into  the  waveguide.  For  waveguides 
with  verv  small  active  regions,  this  poses  a  serious  limitation. 

Despite  these  limitations,  there  is  a  wide  range  of  device  parameters  for 
which  the  technique  is  very  useful.  In  addition  to  the  layered  waveguide  dis¬ 
cussed  here,  the  technique  can  also  be  applied  to  optical  strlpllne  waveguides 
and  stripe  geometry  lasers. 


Lite  of  Illustration* 


Fig.  1.  Experiments!  and  theoretics!  far  field  pattern*  for  a 
multl-laver  laser  structure. 

Fig.  2.  Index  of  refraction  profile  of  the  multi-layer  structure 
used  in  computing  the  theoretical  far  field  patterns  lr. 
Fig.  3.  Thicknesses  of  layers  are  also  shown. 

Fig.  3.  Halfpower  beaavldth  of  fundamental  mode  of  the  structure 
in  Fig.  It  as  a  function  of  An*  for  certain  values  of  An, . 
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field  Solutions  for  the  Lateral  Modes  of 
Stripe  Geometry  Injection  Lasers* 
by 

J.  B.  Delaney  and  J.  K.  Butler 
Electrical  Engineering  Department 
Southern  Methodist  University 
Dallas.  TX  75275 

ABSTRACT 

Tne  mode  properties  of  stripe  geometry  lasers  are  analyzed  from  solutions 
o'  a  bounaerv  technique  which  inextricably  ties  the  vertical  and  lateral  mode 
snaoes.  Assuring  a  parabolic  dielectric  distribution  in  the  active  laver, 
attention  is  'ocused  or  the  lateral  mode  shapes  as  a  function  of  active  layer 
thicknesses;  results  are  compared  to  those  obtained  using  a  simple  Henrite- 
Gaussian  field.  In  particular,  near*fields  of  the  twr  models  tre  compered  and 
di screpancies  in  half-power  *u11  widths  o'  tne  lateral  intensities  are  found 
to  be  greater  than  25*  *or  thin  active  layer  thicknesses.  --0.1  urr,.  The  new 
model,  which  accurately  matches  fields  at  the  vertical  heterolunction  bound¬ 
aries.  predicts  relatively  high  oa'n  values  at  threshold  for  thin  active  lay¬ 
er  structures. 
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Stripe  geometry  AT  Ga^_xAs  lasers  nave  oeer.  modeled  with  a  complex  dielec- 

( 1  -A  1 

trie  constant  in  the  active  layer  of  the  fonrv  ' 

>c(x)  ■  -  k0'ai*?  0) 

where  •:  is  the  complex  dielectric  constant  at  x»C  and  kg  is  the  free  space 
wave  nun&er.  With  this  dielectric  distribution,  modes  can  be  described  by  sim¬ 


ple  Hermite- Gaussian  functions. 


(5-7) 


Such  quantities  as  peak  active  region 


gain  and  lateral  index  changes  can  be  estimatec  from  measurements  of  near/ 
far  field  half-power  widths  and  a  virtual  waist  position.  However,  if  the 
fields  are  matenee  at  the  netercjunctions  (x*^  of  Fig.  1)  It  is  necessary  to 
use  a  summation  o*  Hermite-Gaussian  functions  to  describe  a  mode,  with  the  con- 
seauence  that  the  near  and  ‘tr  radiation  fields  are  altered^  We  briefly 
compare  the  na;or  discrepancies  which  arise  between  using  a  single  term  and  a 
summation  to  describe  the  lateral  'ield. 

Referring  to  the  geometry  of  Fig.  1,  a  summary  of  the  field  representa¬ 
tion  for  a  symmetric  TE  mode  of  the  form  expdwt-yz) ,  where  y»  -  y  4  IB,  is 


Ex(x,y)  •  :  A,cos,cjly)Ht(.e  kQx)exp(-  »  ak^x2) 


!x|<f  (2) 


1 

E^(x,y)  ■  /  B(y)cos(>;x)exp( (|  -  y)(x^  -  k^  y*)^)dx  !x  >w  (3) 

where  (■0,2,4,...  for  even  modes.  H  is  the  Hermite  polynomial  and  r  Is  the 
dielectric  constant  for  the  sandwiching  layers.  The  vertical  eigenvalue,  , 
satisfies 


*0^*0  *  4  kp^a{2l4l) 


(4) 


The  far  field  Intensity  In  the  lateral  plane  is'" 

I  •  1qCos£  |ex  O.kgSine)!2  (5) 

where  e^  Is  the  two  cimensional  Fourier  transform  of  Ex(x,y). 


Comparison  of  the  radiation  field  for  matched  fields  to  the  single  tenr 

case  is  accomplished  with  the  parameter  "a"  In  equation  (1)  defined  as^' 

1  1 

8  ,^‘!(2R+  (6) 

*0*0  u 


6n/6g 


wnere  x*Xq  is  defined  as  the  point  at  which  the  gain  has  fallen  to  zero,  i.e., 
6g»9Q.  6n  is  the  refractive  index  change  at  x»Xq. 

Nume'icel  solutions  illustrating  the  differences  between  the  two  models 
are  presented  in  Figs.  2  and  3.  in  Figs.  2(c)  and  (d)  the  fundamental  mode 
half-power  full  width  B  o*  the  lateral  near- field  is  plotted  against  the  full 
gain  width  2^  for  various  active  region  thicknesses  d.  The  small  cashed  lines 
are  for  the  single  term  case;  the  solid  lines  are  for  the  summation  case  arising 
from  matching  fields  at  the  heterojunctions;  and  the  large  dashed  curve  corres¬ 
ponds  to  the  asymptotic  case  as  d-lerge  (e.g.,  1.0  um).  It  should  be  noted  at 
the  outset  tnat  the  two  models  coincide  for  large  active  regions  and  display 
contrasting  trends  as  d  gets  small  with  less  energy  confined  to  the  active 
layer.  The  beamwidth  Is  a  soft  function  of  d  for  the  summation  case, 
with  8  increasing  as  more  energy  1$  confined  to  the  active  region.  In  contrast, 
beamwidths  for  the  single  term  case  change  more  rapidly  with  B  decreasing 
as  the  active  layer  progresses  to  a  large  thickness.  These  trends  are  graph¬ 
ically  exhibited  in  Fig.  j  (b),  where  B  is  plotted  directly  against  d.  The 


fact  that  the  two  models  coincide  in  the  asymptotic  case  d-large  is  Illustrated 
in  the  expansion  coefficients,  At,  used  in  equation  (2).  As  o-large  anc  all  the 
light  is  confined  to  the  active  layer,  we  find  that  |AJ»1,  lA^.-C  for  LfQ,  and 
the  late-el  field  expression  for  the  sumation  case  is  Gaussian.  As  Fig.  3(b) 
indicates,  the  difference  in  the  two  models  is  negligible  for  d>0.5jir.. 

The  far-field  half  power  full  widths, f  »  are  calculated  for  completeness 

if 

and  presented  in  Figs.  2  (a)  and  (b)  and  Fig.  3  (a).  Again,  the  two  models 
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overlay  for  d>0. Sum.  A  value  of  R«-1C  includes  a  small  refractive  Index  rise 

at  x«x^  over  the  value  at  x«0.  For  instance,  with  2x^»12um,  d«D.2um,  £n«0.2. 

the  peak  gain  numerically  calculated  for  the  summation  case  is  Sq"135  cm‘\ 

This  is  the  gain  required  to  propagate  a  mode  with  a  modal  gain  G^SOcm'1  to 

satisfy  radiation  end  losses.  Thus  6n  •  4gR  •  -0.000135. 

•  • 

Typical  peak  gains  are  plotted  against  the  active  region  thickness  in  Fig .4 
for  the  case  R«0,2xr»  12um.  The  solid  lines  are  for  the  sumation  and  the  dash¬ 
ed  lines  are  for  the  single  term..  The  differences  in  peak  gain  values  are  quite 
notable  for  d<0.2ym.  The  peak  gain  follows  similar  trends  for  other  values  of 
the  parameter  R.  The  more  rapid  rise  of  the  gain  as  d  drops  in  the  solid  curves 
implies  higher  currents  at  threshold. 

A  worst  case  example  of  the  difference  in  the  two  models  is  R«-10*^, 
2xQ»12um,  d»0.1um.  in«0.2.  The  single  term  near-field  spot  width  is  29t  larger 
than  that  for  the  summation  model.  The  single  term  far-field  width  is  27i 
smaller  than  its  summation  counterpart. 

In  conclusion,  a  simple  Hermite  Gaussian  mode  description  of  the  lateral 
field  has  the  inherent  attractiveness  of  analytic  simplicity,  with  the  concomi¬ 
tant  ease  of  matching  experimental  evidence  of  the  laser.  However,  the  approxi¬ 
mation  loses  accuracy  for  very  thin  active  layers.  Properly  matched  fields 
yield  an  expression  for  the  lateral  field  intensity  which  can  significantly 
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differ  froir  a  single  term  description  as  the  active  layer  narrows.  Resorting 
to  the  more  complicated  sumnation  expression  has  the  disadvantage  in  that  a 
computer  has  to  be  used  to  manaoe  an  Infinite  series.  When  the  active  layer 
is  thicker  than  0.25uir.,  It  Is  likely  that  the  error  In  using  one  term  for  the 
lateral  field  is  not  more  crucial  than  the  error  unavoidable  in  laboratory  ex¬ 
periment. 
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FIGURE  CAPTIONS 


Figure  1.  Cross-section  of  the  active  region  of  a  stripe  geometry  laser.  The 
complex  dielectric  constant  of  the  active  layer  is  parabolic  in  x. 

Figure  2.  B  Is  the  fundamental  mode  half -power  full  width  of  the  lateral  near¬ 
field  Intensity.  Op  is  the  fundamental  mode  half-power  full  width  of  the  lateral 
far  field  intensity.  2x^  is  the  full  gain  width  of  the  active  layer  under  the 
stripe.  Small  dashed  lines  correspond  to  a  single  term  representation  of  the 
lateral  field;  solid  lines  are  for  the  lateral  field  represented  by  a  simulation 
of  Henri te-Gaussian  functions;  and  the  large  dashed  lines  are  for  the  asymptotic 
case  d-large  (1.0  ;m>). 

Figure  3.  B  and  e,,a re  plotted  against  the  active  region  thickness  d  in  microns. 
The  dashed  lines  are  for  the  case  where  the  lateral  field  is  a  single  term. 

The  solid  lines  are  for  a  simulation  of  Hermite-Gaussian  functions.  2xg»12tfn. 

All  curves  are  'or  the  fundamental  mode. 

Figure  A.  Peak  gain  is  plotted  against  the  active  region  thickness  d  In 
microns  for  the  fundamental  and  second  mode.  The  solid  curves  are  for  a  sum¬ 
mation  while  the  dashed  curves  are  for  a  single  term  representation  of  lateral 
mode. 
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ELECTROMAGNETIC  MOOES  Of  HETEROJUNCTION  INJECTION  LASERS 
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This  research  la  concerned  with  the  theoretical  nodding  of 
he teroj unction  Injection  lasers.  Studies  have  been  aade  on  two 
dimensional  radiation  patterns,  mode  reflectivities  due  to  cleved 
facets,  transverse  node  selection,  and  the  relationships  between  the 
threshold  gain  and  aluminum  concentration  of  various  grown  layars  in 
the  CaAs-AlCaAs  device  structure. 

e 

Two  dimensional  mods  confinement  Is  discussed  for  several  laser 
structures  Including  stripe  geometry  devices.  Radiation  patterns  and 
mode  ref lectivltlea  were  calculated  using  the  two  dimensional  model. 

Transverse  mode  selection  due  to  a  thin  film  costing  on  the 
laser  facet  is  discussed. 

The  threshold  behavior  for  C W  operation  Is  investigated.  The 
aluminum  mola  fraction  In  the  active  region  determining  the  lasing 
-avelangth,  affects  the  active  region  gain  at  threshold.  Lasers  with 
a  given  active  region  width  and  Index  step  have  the  smallest  threshold 
gain  when  aluminum  fraction  of  the  active  region  Is  the  largest;  the 
peak  -fractional  value  of  aluminum  cannot  exceed  0.37. 
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SATIATION  FIELDS  OF  DIELECTRIC  OPTICAL  WAVEGUIDES 

Advisor,  Professor  Jerome  K.  Butler 
Doctor  of  Philosophy  degree  conferred:  August  6,  1979 
Dissertation  completed:  July  3,  1979 

School  of  Engineering  and  Applied  Science,  Southern  Methodist  University 

At.  experimental  method  of  measuring  the  far  field  radiation 
partem  cf  dielectric  optical  waveguide*  i*  presented.  The  patterns 
of  three-layer  and  five-layer  structures,  obtained  using  this  technlcue, 
are  given.  Matching  the  experimental  results  with  patterns  calculated 
using  the  theory  cf  layered  waveguides  determines  the  geometry  and  index 
of  refraction  profile  of  the  waveguide. 

Dispersion  characteristics  and  radiation  fields  cf  an  optical 
stripline  waveguide  radiating  into  free  space  are  *alcuiated.  The 
naif power  beaswicch  cf  the  fundamental  mode  in  a  plane  parallel  to  the 
dielectric  layers  is  used  as  an  indication  cf  the  degree  of  optical 
confinement  under  the  stripe.  Flots  of  the  beaou'idth  versus  the 
various  waveguide  parameters  indicate  the  parameter  values  which  give 
the  best  optical  confinement.  An  approximate,  and  computationally 
simpler ,  method  of  calculating  this  beaawidth  la  also  discussed. 
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Injection  Laser 
Radiation  Patterns 
Waveguide  Modes 
Mode  Stability 


This  research  was  concerned  with  (1)  the  electromagnetic  modes  and  (2>  the 
radiation  pattern  characteristics  of  contesmorarv  injection  laser  structures. 
The  objective  was  to  Investigate  laser  structures  which  can  produce  stable- 
beaa  radiation  patterns.  The  control  of  lateral  fields  in  lasers  has  been 
accomplished  with  manv  novel  devices  such  as  the  channeled-subatrate-planar, 
burled  heterostructure  and  strip-loaded  waveguide.  The  drawback  to  the  atrlpe 
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laser  Is  chat  It  Is  laterally  unstable.  Much  of  our  work  has  been  alaed  at 
rigorously  modeling  the  electromagnetic  modes  of  stripe  lasers.  We  have  also 
applied  our  analysis  to  other  lasers  such  as  the  channeled-substrate-planar 
device. 


